ALGEBRAS, GROUPS AND GEOMETRIES Volume 40, Number 1, June 2024 


GROUPS 
AND 


GEOMETRIES 


FOUNDED IN 1984. SOME OF THE PAST EDITORS 
INCLUDE PROFESSORS: A. A. ANDERSEN, 
V. V. DEODHAR, J. R. FAULKNER, 

M. GOLDBERG, N. IWAHORI, M. KOECHER, 
M. MARCUS, K. M. MCCRIMMON, R. V. MOODY, 
R. H. OHEMKE, J. M. OSBORN, J. PATERA, 
H. RUND, T. A. SPRINGER, G. S. TSAGAS. 


EDITORIAL BOARD 


M.R. ADHIKARI 
A.K. ARINGAZIN 

P. K. BANERJI 

A. BAYOUMI 
QING-MING CHENG 
A. HARMANCY 

L. PR HORWITZ 

C. P. JOHNSON 

S.L. KALLA 

N. KAMIYA 

J. LISSNER 

J. LOHMUS* 

R. MIRON 

M.R. MOLAEI 

P. NOWOSAD 
KAR-PING SHUM 
SERGEI SILVESTROV 
H. M. SRIVASTAVA 
G.T. TSAGAS* 


* In memoriam 


FOUNDER and 
Editor In Chief 
R. M. SANTILLI 


ALGEBRAS, GROUPS AND GEOMETRIES 


Established in 1978 by Prof. R. M. Santilli at Harvard 
University Hadronic Journal and Algebras, Groups 
and Geometries have been regularly published since 
1978 without publication charges and are among the 
few remaining independent refereed journals. 


This Journal publishes advances research papers 
and Ph. D. theses in any field of mathematics 
without publication charges. 


For subscription, format and any other information 
please visit the website 
http://www.hadronicpress.com 


HADRONIC PRESS INC. 
35246 U. S. 19 North Suite 215 
Palm Harbor, FL 34684, U.S.A. 
http://www.hadronicpress.com 

Email: info@hadronicpress.com 
Phone: +1-727-946-0427 


ALGEBRAS, GROUPS AND GEOMETRIES 


Founder and Editor in Chief 
RUGGERO MARIA SANTILLI 


The Institute for Basic Research 
35246 U. S. 19 North Suite 215, Palm Harbor, FL 34684, U.S.A. 
Email: resarch@i-b-r.org; TEL: +1-727-688-3992 


R.M. ADHIKARI Calcutta Mathematical 
Society AE-374, Sector-1, Salt Lake City 
Kolkata, West Bengal 700064 

A.K. ARINGAZIN Department of 
Theoretical Physics, Institute for Basic 
Research, Eurasian National University, 
Astana 010008, Kazakhstan 
aringazin@mail.kz 


P.K. BANERJI Department of Theoretical 


Physics, Faculty of Science, I.N.V. 
University, Jofhpur-342 2005, India 
bannerjipk@yahoo.com 

A. BAYOUMI Alazhar UniversityFaculty 
of Science Mathematics Department, 
Naser City 11884 CairoEgypt 
draboubakr@yahoo.com 

QING-MING CHENG Department of 
Mathematics, Faculty of Science Josat 
University, Sakado, Saitama 350-02 Japan 
cheng@math.josai.ac.jpbr 

A. HARMANCY Department of 
Mathematics, University of Hacettepe, 
Beytepe Campus, Ankara, Turkey 
harmaci(@hacettepe.edu.tr 

L. P. HORWITZ Department of Physics, 
Tel Aviv University, Ramat Aviv, Israel 
horwitz@taunivm.tau.ac.il 

C. P., JOHNSON Department of 
Mathematics, Mississippi State University 
P. O. Box MA, Mississippi State, MS 39762 
cjohnson@math.msstate.edu 

S. L. KALLA Department of 
Mathematics, Vyas Institutes of 

Higher Education, Jodhpur 

342008 India 

shyamkalla@gmail.com 

N. KAMIYA University of Aizu 

Center for Mathematical Sciences, 
Tsuruga, [kki-machi Aizu- 

Wakamatsu City Fukushima, 965- 

8580 Japan 


J. LISSNER, Alumnus 
Foukzon Laboratory 

Center for Mathematical Sciences 
Technion-Israel Institute of Technology 
Haifa, 3200003, Israel 

L. LOHMUS* Estonia Academy of 
Sciences, 142 Rita Street, Tartu 
202400, Estonia 

R. MIRON Faculty of Mathematics 
University “Al. I. Cuza” 6600 Isai, 
Romania rmiron@uaic.to 

M.R. MOLAEI Department of 
Mathematics, Universty of of Kerman, 
P.O. Box 76135-133 Kerman, Iran 
mrmolaei @maill.uk.ac.ir 

P. NOWOSAD Mathematics 
FFCLRP - USP,Univ. of San 

Paulo 14040-901 Ribeirao Preto 

SP, Brazil 

KAR-PING SHUM Institute of 
Mathematics, Yunnan University, 
Kunming, China 
kpshum@ynu.edu.cn 

S. SILVESTROV School of 
Education, Culture and 
Communication (UKK) Malardalen 
University, Box 883, 71610 

Vasteras, Sweden 
sergei.silvestrov@madh.se 

H. M. SRIVASTAVA Department of 
Mathematics and Statistics, 

University of Victoria, Victoria, B.C. 
V8W 3P4, Canada 
hmsri@uvvm.uvic.ca 

G. EF TSAGAS* Mathematics 
Department, Aristotle University 
Thessaloniki 54006, Greece 


*In Memorian 


yZ0Z euns ‘| JOqUINN ‘OV eWN/OA SAINLAWOAD GNV SdNOUD ‘SvydsN1V 


INC. 


Vv 


vi 
Ww 
lw 
oc 
a. 
~ 
2 
O 
cc 
Q 
< 
<= 


0741 - 9937 


few remaining independent refereed journals. 
This Journal publishes advances research papers 
and Ph. D. theses in any field of mathematics 
without publication charges. 
please visit the website 
http://www.hadronicpress.com 
HADRONIC PRESS INC. 
35246 U. S. 19 North Suite 215 
Palm Harbor, FL 34684, U.S.A. 
http://www.hadronicpress.com 
Email: info@hadronicpress.com 
Phone: +1-727-946-0427 
ISSN 


University Hadronic Journal and Algebras, Groups 
and Geometries have been regularly published since 


Established in 1978 by Prof. R. M. Santilli at Harvard 
1978 without publication charges and are among the 
For subscription, format and any other information 


ALGEBRAS, GROUPS 
AND GEOMETRIES 


VOLUME 40, NUMBER 1, JUNE 2024 


OPERATOR SEMIGROUPS ACTING ON A I-SEMIGROUPS, 1 
Johnson Aderemi Awolola 
Department of Mathematics 
Federal University of Agriculture 
Makurdi, Nigeria 
Musa Adeku Ibrahim 
Department of Mathematics 
Federal University Lokoja 
Kogi State, Nigeria 
AN APPLICATION OF SMITH NORMAL FORM, 11 
Haohao Wang 
Department of Mathematics 
Southeast Missouri State University 
Cape Girardeau, MO 63701 
ABSTRACT PETER-WEYL THEORY FOR SEMICOMPLETE 
ORTHONORMAL SETS, 23 
Olufemi O. Oyadare 
Department of Mathematics 
Obafemi Awolowo University 
Ile-Ife, 220005, Nigeria 
SERIES ANALYSIS AND SCHWARTZ ALGEBRAS OF SPHERICAL 
CONVOLUTIONS ON SEMISIMPLE LIE GROUPS, 41 
Olufemi O. Oyadare 
Department of Mathematics 
Obafemi Awolowo University 
Ile-Ife, 220005, Nigeria 
QUANTUM RELATIVISTIC EQUATION FOR A 
PROBABILTY AMPLITUDE, 61 
Yu.M. Poluektov 
National Science Center 
“Kharkov Institute of Physics and Technology” 
61108 Kharkov, Ukraine 


SINGLE TERM APPROXIMATIONS OF HYPERGEOMETRIC FUNCTIONS, 79 

Leda Galue 

Centro de Investigacion de Matematica Aplicada 

Fac. De Ingenieria, Universidad del Zulia 

Apartado 10.182, Maracaibo, Venezuela 

Shyam L. Kalla 

Kuwait University, Fac. of Science, Dept. of Mathmatics 

P.O. Box 5969, Safat 13060, Kuwait 

Virginia Kiryakova 

Institute of Mathematics, Bulgarian Academy of Sciences 

P.O. Box 373, 1090 Sofia, Bulgaria 
FRACTALS GENERATED BY MOBIUS TRANSFORMATIONS AND SOME 
APPLICATIONS, 103 

Sergiu Corbu and Mihai Postolache 

“Politehnica” University of Bucharest 

Department of Mathematics I 

Splaiul Independentei 313 

77206 Bucharest, Romania 
A SERIES-ITERATION METHOD IN THE THEORY OF ORDINARY 
DIFFERENTIAL EQUATIONS, 113 

D. Dimitrovski 

Institute of Mathematics 

Faculty of Science, Skopje 

M. Mijatovic 

Institute of Physics, Faculty of Science, Skopje and 

Institute of Nuclear Sciences 

Vinca, Belgrade, Yugoslavia 


ALGEBRAS, GROUPS AND GEOMETRIES 40 1-10 (2024) -l- 
DOI 10.29083/AGG.40.01.2024/SC1 


OPERATOR SEMIGROUPS ACTING ON A I'-SEMIGROUPS 


Johnson Aderemi Awolola 
Department of Mathematics 
Federal University of Agriculture 
Makurdi, Nigeria 
remsonay@yahoo.com 


Musa Adeku Ibrahim 
Department of Mathematics 
Federal University Lokoja 
Kogi State, Nigeria 
musa.ibrahim@fulokoja.edu.ng 


Received August 15, 2023 


Abstract 
The concept of I-semigroup is a generalization of semigroups. In this pa- 
per, we briefly introduce the action of left (right) operator semigroups on 
a I-semigroup and deduce in particular that there exists an inclusion pre- 
serving bijection between the set of all right ideals of S and the set all right 
ideals of LD x S. 
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1 Introduction 


The notion of I in algebra was first introduced by Nobusawa [8] as a gen- 
eralization of ring in the form of I-ring. Let M and I be additive groups 
such that for all a,b,c € M and 7, 8,a € T, we have ayb € M and ya8 €T 
for every a,b, y and G, then M is called a I-ring if the following conditions 
are satisfied: 


(i) (a1 + a2)yb = ayyb + aeyb, 
aly +2)b = ayib + ay2b, 
_ay(b; + bz) = ayb; + abn, 


(ii) (a7b)Be = ay(bBc) = a(yb8)c, 
(iii) if ayb = 0 for any a,b € M, then y = 0. 


The structure of I'-rings as initiated by [8] was studied by Barnes [{1], Luh 
(7], Ravishankar and Shukla [9], Buys and Groenewald [3], Booth [2] and 
Kyuno [6]. Motivated by this generalization of a ring, Sen [13] defined the 
concept of [-semigroup. Later, Sen and Saha [14] redefined the [’-semigroup 
by weakening slightly the defining conditions of [-semigroup to ensure it 
preserves semigroup structure. The development of I-semigroups hinges on 
the fact that subsets of a semigroup naturally inherits associativity but not 
necessarily closed. As a result of this, various generalizations and analogues 
of corresponding results in semigroup theory have been obtained based on 
the modified definition (see [10, 15, 16, 17]). 

In an attempt to broaden the theoretical aspect of [-semigroup the- 
ory, Dutta and Adhikari [4] slightly changed the defining conditions of I’- 
semigroup by Sen and Saha [14] and then introduced the notion of left 
operator semigroup and right operator semigroup of a I’-semigroup. 

In [4], the authors described the relationship between the set of I-ideals 
and operator semigroups. In relation to the concept, Dutta and Chattopad- 
hyay [5] initiated the notions of uniformly strongly prime semigroup, uni- 
formly strongly prime ideal, Rees congruence on I’-semigroup and uniformly 
strongly prime ['-semigroup and studied these via operator semigroups of 
[-semigroup. Sardar et al. [12] showed that the left operator and right 


operator semigroups of a I’-semigroup with unities are Morita equivalence 
monoid and further established that there is a close connection between the 
Morita equivalence of monoids and I’-semigroups. Although there is signif- 
icantly number of published results in literature on operator semigroups of 
a I-semigroup, however the aspect of operator semigroups acting on a [- 
semigroup observed by [11] has not been given much attention. This serves 
as a motivation to write this paper and we deduce some results. 


2 Preliminaries 


We recall some definitions and results related to this paper. 


Definition 2.1 Let S and T be two non-empty sets. Then S' is called a 
['-semigroup if there exist mappings S xT x S + S$ | (a,a,b) -+ aabe S 
andT x SxT —T | (a,a,8) --+ aaG € I satisfying the identities 
aa(bBc) = a(abB)c = (aab) Bc for all a,b,c € S anda, BET. 


The modified definition of !-semigroup by Sen and Saha [14] may be 
regarded as one-sided I’-semigroup. 


Definition 2.2 Let S = {a,b,c...} andl = {a,6,7,...} be two non-empty 
sets. Then S is called al'-semigroup if there exists a mapping S xT x S —> 
S| (a,a,b) —> aab € S satisfying the property (aab)Bc = aa(bBc) for all 
a,b,cE S anda, Be). 


A [-semigroup S$ is called commutative if cay = yax for every x,y € S 
andaeél, 
Let A and B be two subsets of a I-semigroup S. We define the set 


ADB = {ayb|aeé A, be Band yeET}. 


For simplicity we write a'B, AI'b and AyB instead of {a}I'B, AI'{b} 
and A{y}B respectively. 

Let S be an arbitrary semigroup and I’ be a non-empty set. Define a 
mapping S x I x S —> S by aab = ab for all a,b € S anda eT. It is easy 
to see that S is a I-semigroup. Thus, a semigroup can be considered as a 
I’-semigroup. 

In the following, some examples of I'-semigroups are presented. 


Example 2.1 Let S = Z be the set of all integers andl = {n|ne€ N}. 
Then S is a T'-semigroup with the operation defined by aab =a+a+b for 
alla,be S anda él. 


Example 2.2 Let S be a set of all negative rational numbers. Obviously S 


is not a semigroup under usual product of rational numbers. Let T = {—+ 


P 
p is prime}. Let a,b,c € S anda,B €T. Now if aab is equal to the usual 
product of rational numbers a,a:,b; then aab € S and (aab)Bc = aa(bBc). 


Hence S is a V-semigroup. 


Example 2.3 Let S = {—i,1,0} andl =S. ThenS is al-semigroup with 
respect to multiplication of complex numbers whereas S does not reduce to 
semigroup with respect to multiplication of complex numbers. 


The following definitions and theorems can be found in [5] except Defi- 
nitions 2.3. 


Definition 2.3 Let S be al'-semigroup. A nonempty subset A of S is called 
left (right) T'-ideal of S if STA C A (ATS C A). Further, a non-empty A 
of aI'-semigroup S' is called I’-ideal if A is both a left and a right T’-ideal of 
S. 


Definition 2.4 Let S be al'-semigroup. Let L and R be the left and right 
operator semigroups of the [-semigroup S. If there exist an element [e, 5] € 
L ((6,e] € R) such that eér = x (xbe = 2) for allx € S, then S is said to 
have the left unity [e, 6] (right unity [6, el). 


Definition 2.5 Let S be a T'-semigroup. We define a relation p on S xT 
as follows: 


(x, a)p(y, B) <> ras = yBs, Vs ES. 


Obviously p is an equivalence relation. Let [x,a] denote the equivalence 
class containing (a,a). Let L = {[z,a]): 2 €S, a€T}. Then L is a 
semigroup with respect to multiplication defined by {x, al[y, B] = [xay, ]. 
The semigroup L is called the left operator semigroup of S. Similarly, the 
right operator semigroup R of aT-semigroup S is defined as R = {[a, a] : 
ae, c€ S}, where [a,z][6,y] = [a,rGy], for allz,y eS anda, BET. 


Example 2.4 Using Example 2.3, 

SxP= {(—4, i), (0, 0), (4, 2), (—i, 0), (—1,2), (0, —4), (0,7), (4, —t), (4,0) } 
and let p = {(—1, i), (0,0), (i, 2), (—é, 2), (¢, -2)} be @ relation on S xT. 
By routine calculation, it is obvious that p is an equivalence relation. The 

equivalence class [x, a] = {(y, 8) € p | (y, B)p(z, a)}. Therefore, 
[-4, —1] = {(—4, —i), (4, 4)} 
* (0,0) = {(0,0)} 
[é,2] = {(i, ‘); (—é,-7)} 
[—4,¢] = {(—#, ‘); (3, —A)} 
[2, —1] = {(, —t), (2, =4)} 
=> [-i, —7] = [2,7] and [-i, i] = [i, —é] and the set form L = {[—4, —4], [2, a], [-7, 2], [2, —2]} 
is a left operator semigroup of S. Similarly, the right operator semigroup R 
of S can be obtained. 


Theorem 2.1 Let S be aI’-semigroup. If le, 6] is left unity of S, then [e, 6] 
is the identity element of L. 


Theorem 2.2 Let S be aT-semigroup. If |u, f] is right unity of S, then 
(ut, f] 1s the identity element of R. 


Theorem 2.3 Let S be aIT-semigroup with left and right unities and L be 
its left operator semigroup. 


(i) If Q is aT-tdeal of S', then gt is a T'-ideal of L. 

(ii) If P is aT -ideal of L, then P* is aT-ideal of S. 
Theorem 2.4 Let S be al'-semigroup with left and right unities and R be 
its right operator semigroup. 

(i) IfQ is aT -ideal of S, then gt is a T’-ideal of R. 

(ii) If P is aT -ideal of R, then P* is a T-ideal of S. 


Theorem 2.5 Let S be a T'-semigroup with left and right unities and let 
L and R be its left operator semigroup and right operator semigroup re- 
spectively. Then there is an inclusion preserving bijection between the set 
of all T’-ideals of S and the set of all T'-ideals of L(R) via the mapping 


Oot (Q—+ Q*) where Q is a T-ideal of S. 


3 Main Results 


Throughout S stands for one-sided I’-semigroup unless otherwise mentioned. 


The following proposition and theorem show the commutativity and 
isomorphism of operator semigroups. 


Proposition 3.1 Let S be a commutative T-semigroup anda € TV. Then 
the left operator semigroup L and the right operator semigroup R of S are 
commutative. 


Proof 
The proof is straightforward. 


Theorem 3.1 [f S is a commutative T'-semigroup and a € I, then the op- 
erator semigroup L and the right operator semigroup R of S are isomorphic. 


Proof 
Define f : L — R by f([a,a]) =[a,a]. Let [a,a] = [b,a]. Then aas = bas 
for all s € S. Since S is commutative, saa = sab for all s € S. So, 
(a, a] = [a,b]. Thus, f is well-defined. The mapping is injective, since 
[a,a]= [a,b] => saa=sab VsES 
=> aas= bas 
=> [a,a] = [b,a] 
Again, f is surjective, since for any [a,a] € R,a € S anda ET, we 


have [a, a] = f[a,q]. 
Also, for alla,be S andaeT, 


Ff ([a, a][b, a]) = f([aab,a]) = [a, aad] 


I 


[a, baa] 
[a,b][a,a] 

[a, a][a, 5] 
F(la,a]) f({oe]) 


So, f is ahomomorphism. Therefore, L and R are isomorphic. 


I 


In the following we consider operator semigroups acting on a I’-semigroup. 


Definition 3.1 Let S be al'-semigroup and its left and right operator semi- 
group are respectively L and R. Then for everya€e S anda ET, we define 
LxS—S andS x R—} S respectively as follow: 


(2, a]s := aas and s[a, a] := saa. 
The following remark follows from Definition 3.1. 
Remark 3.1 Jf S' is a commutative [-semigroup, then Lx S=S x R. 


Example 3.1 Clearly, Example 2.3 shows that S' is commutative and from 
Example 2.4 it is easy to verify that Lx S = S x R for every s € S and 
aeéelr. 


Proposition 3.2 If A is an ideal of S, then L x A is an ideal of L x S. 


Proof 

Suppose that A is an ideal of and [a, alt € L x A. Then for every s,t € S 
and 8 € T, we have aa(sft) € A. Thus, [a,a][s,8]t = [aas, Blt € Lx A. 
Similarly, we can prove that [s, 6][a, o]t = [s8a,a]t € L x A. Hence, L x A 
is an ideal of L x S. 


Theorem 3.2 There exists an inclusion preserving bijection between the 
set of all right ideals of S and the set all right ideals of L x S. 


Proof 

Suppose that Z(S') and Z(L x) sre the sets of all right ideals of S and Lx S 
respectively. Clearly, the mapping f : Z(S) —> Z(L x S) by f(I) =LxI 
is well-defined. Since J is a right ideal of S, JTS C J. Thus, CL x I. 
On the other hand, since S has a left unity, Lx J C(Lx IT'S CI. Thus, 
LxI=TI. Hence, f is bijective. It remains to show that f is an inclusion 
preserving mapping. Let J and J be two right ideals of S such that J C J. 
We have to show that L x I CL x J. Let [a,a] € L. Then for every s € S 
we have s € J and so [a,a]s € L x I. Since J C J, we have s € J. Hence, 
[a,a]s € L x J. Therefore, L x I C Lx J. Hence the proof. 


Remark 3.2 Similar characterisations can be proved for right operator 
semigroup acting on a T'-semigroup and right T’-ideal. 
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Abstract 


Modules are a generalization of the vector spaces of linear alge- 
bra over a ring instead of over a field of scalars. Some of the re- 
sults in linear algebra are extremely useful in studying modules. In 
this paper, we highlight the use of Smith normal form in studying 
finitely generated modules over a univariate polynomial ring. As an 
application, we take advantage of Smith normal form to understand 
the structure of certain syzygy modules. 
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1 Introduction 


Linear algebra deals with vector spaces and linear transformations over 
a field such as real or complex numbers. Linear algebra is very well un- 
derstood, serves as a fundamental tool for mathematical economics, data 
science, machine learning, and has many applications, from mathemati- 
cal physics to modern algebra and coding theory. 

Modules are a generalization of the vector spaces of linear algebra 
over a ring instead of over a field of scalars. A fundamental fact of lin- 
ear algebra over a field is a finitely generated vector space has a basis — 
the minimal generating (spanning) set of a vector space are linearly inde- 
pendent and therefore form a basis. However, modules are more compli- 
cated than vector spaces; for instance, not all modules have a basis. Only 
a special family of modules called free modules have a basis. It remains 
an interesting and an active research area to find the minimal generating 
set, or an upper bound for the minimal generating set for different kinds 
of modules under various of conditions [2], [3], [4], and [6]. 

In linear algebra, matrix factorization is a mathematical technique 
used to decompose a matrix into the product of multiple matrices. There 
are various ways to decompose a matrix depending on the context in 
which that matrix is used. For instance, in computer science and machine 
learning, singular value decomposition is one of the most important com- 
putational method. Whereas, representing a matrix in the Smith normal 
form has many computational applications in number theory, group the- 
ory, and homological algebra. 

In this short paper, we focus on the Smith normal form, and its ap- 
plications in modules over univariate polynomial rings K[z|] with K an 
algebraically closed field of characteristic zero. Our goal is to find a min- 
imal set of generators for certain syzygy modules using Smith normal 
form. 

Recall that given a generating set {f1,..., fm} of an ideal (or a module) 
over a ring R, a relation or the first syzygy between the generators is an 
m-tuple (a1,...,@m) € R™ such that 


ay fy, + ao fo +--+ Omim = 9. 


ae 


The set of syzygies form a syzygy module, denoted by S = Syz(fi,...5 fi im) 
A Koszul syzygy is one of the form 


(fifi + (—fi) fi = 0, fori A 3. 


The Koszul syzygies generate a submodule of the syzygy module S. 
Our attention is centered on the submodule K that are generated by 
the “Koszul-like” forms over a univariate polynomial ring K[2]: 


ie ex ,0 a L<i<j<m. 
gcd( fi, f5) = ~ 


a 


O.eas Oise 


0 eee are 
ged By 
We seek the answers to the question that how the minimal sets of gener- 
ators for K and S are related. 

This paper is structured as the following. We begin in Section 2 with a 
brief review of results concerning the Smith normal form over Principal 
Ideal Domains (PIDs), and the Hilbert-Burch Theorem. We then study 
the minimal set of generators for the submodule K of a syzygy module 
S in Section 3. The main results of this paper are Theorems 3.1 and 3.2, 
where we provide an explicit equation to describe a relationship between 
a basis for K and a basis for S. We flush out our theorems by a simple 
and illustrative example. 


2 <A Brief Review 


The Smith Normal Form (SNF) is a canonical form to which a matrix 
can be transformed using elementary row and column operations. It is 
particularly useful when dealing with matrices whose entries come from 
a ring with special properties, such as PIDs. Theorem 2.1 is a well-known 
results concerning SNF over PID. 


Theorem 2.1. (Smith Normal Form over PID [1, Theorem 3.1]) Let R 
be a PID and let A € Mmm(R). Then there is aU € GL(m,R) and a 
V €GL(n, R) such that 


Dy 0 
vay =|" ile 


olde 


where r = rank(A), and D is a diagonal matriz with diagonal entries 
$1, S9,...,8, With s, 40 fori=1,...,7, and 5 | 41 fori=1,...,r—L. 
Additionally, s;’s are unique up to multiplication by a unit and are called 
the elementary divisors, invariants, or invariant factors, which can be com- 
puted (up to multiplication by a unit) as s; = ety where d;(A) is called 
i-th determinant divisor that equals the greatest common divisor of the de- 
terminants of alli xi minors of the matriz A and do(A) := 1. Furthermore, 
if R is a Euclidean domain, the matrices U,V can be taken to be a product 
of elementary matrices. 


In the context of modules (generalizations of vector spaces) over a 
PID, the SNF reveals important information about the structure of the 
module and its relationship to other modules. Throughout this paper, 
we should use Theorem 2.1 to investigate the minimal set of generators 
for modules. To our advantage, we refer SNF as the matrix factorization 
that decomposes A into the product of invertible matrices with a diagonal 
matrix, and write A = U eo 0 

For the convenience of our readers, we also cite a part of the Hilbert-— 
Burch Theorem below. The complete statement of the Hilbert-Burch The- 
orem is beyond the scope of this paper, we only cite the portion of the 
theorem that concerns the first syzygy module Syz(fi,..., fm) over K[z]. 


Theorem 2.2. (Hilbert-Burch Theorem [5]) If M is a module minimal 
generated by m elements over a polynomial ring K(x] over a field K, then 
the first syzygy module of M is always a free module of rank m — 1. 


4 V for some invertible matrices U,V. 


We want to emphasize that this result states that the first syzygy mod- 
ule Syz(fi,--., fm) has a basis consisting of m — 1 linearly independent 
generators over univariate polynomial rings; however, it is not true in 
general for multivariate polynomial rings. 


3 Main Results 


Theorem 3.1. Let {p1,...,Dm—1} be a basis for the syzygy module S = 
Syz(fi,..-)fm) where fi,...,fm €. B= K|z]. Let K be the submodule 


-15- 
generated by the syzygies of the form 


Si fi } Or 
Or a, 0,093 are eG »lsi<jsm 
ged(fi, fi) god( fi, fy) : 


Then 


1. K is minimally generated by m — 1 syzygies of such form, that is, 
rank(K) =m-—1. 


2. There exist at most m —1 distinct Ci,...,Cm—1 € R where C; is the 
smallest degree element in R such that C;p; € K. 


8. The colon ideal (K : S\ = {r € R| rS C K} is the principle ideal 
generated by C = LCM(C),...,Cm-1), the least common multiple of 
C1, - ‘ vy Oey that 18, (K : S) = (LCM(Ci, wae Cat) = (C), 


4. S/K & R/(ti) x R/(te) x +++ x R/(tm—1), where ti,...,tm—1 are the 
elementary divisors of T for some T such that K = ST. 


Proof. Let K be asyzygy submodule generated by the columns of the matrix 


ra fm —Im- my dries os Bus = 
gcd( fi fm) Sata ) 0 fs 0 
~jim T~JIm-1 
: 0 ae ee CD 
0 0 0 0 sani 0 


epee ee et eel 

0 ged(fi,fm—1) 0 gcd(fe,fm—1) gcd(fimsfm—1 
satiny 0 iss ati 0 i Frnt 
gcd(fi,fm) gcd( fo, fim gcd(fmfm—1). 


This matrix is of size m x 1 where | = (3). To simplify the the matrix 
expression, let K; for 7 =1,...,1 be the i-th column of the above matrix. 

Proof of item 1. It is easy to observe that the (m — 1) x (m — 1) 
submatrix formed by the first m — 1 columns and the last m — 1 rows has 
a non-zero determinant. Hence the rank of this matrix is at least m — 1, 
ie., rank(K) > m-—1. On the other hand, since K is a submodule of 
S, rank(K) < rank(S), and rank(S) = m—1 by Theorem 2.2. Thus, 
rank(K) =m-—1. 
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Proof of item 2. Without loss of generality, we may select m — 1 
columns of the matrix that generate the submodule K, and name these 
generators Ky, Ko,...,Km—1. Let a = [a1,...,@m]” € S. Since the matrix 
(Ki, Ko,...,Km—1] € Mmm-i(R) is of rank m — 1, the matrix equation 
[Ki,..., Km-1]X = a has a unique solution over the function field of R, i.e. 


X = [bi /c1, bs/ca, as sbacct/ee4|" € (K(x), with gcd(b,, ¢) =1. 


Let ¢c = LCM(c1,¢2,...,Cm—1) € K[g], set cj = 2 for? =1,...,m—1. Then 
c, € K[x], cX = [eybi, hbo, ..., 0, yom}? € (K[a])™, and 


ca = e([K4, Seay Kyn-1|X) _ (Ki, soe  Km-1](cX) 
C4 by 
CAD) linet F 
= [K},Ko,...,Km-s| = So (eb) Ki eK, 
i=1 
Cip—10m~t1. 
Note since c = LCM(cy,ce,...,€m—1), up to a constant multiple, c is the 


smallest degree element in R such that ca € K. That is, for any factor 6 of 
c that is not a constant multiple of c, Ba ¢ K. 

Since {p1,...,Dm_i} is a basis for S, by the above argument, there exists 
at most m—1 distinct C,, C2,...,Cm—1 € R where C; is the smallest degree 
element in R such that Cyp; € K for each 1. 

Proof of item 3. Let C = LCM (Cj,...,Cm-_1), then Cp; € K for each 
i yields that Cp; € K. Thus, for any a € S, a = $ipi +++ Bm-1pm—1 for 
some (1,...,8m—-1 € R, and 


Ca = C(61P1 + +++ + Bm—1Pm-1) = Bi(C pi) + +++ + Bm-1(CPm-1) € K. 


Therefore, we have (K : S) = (C) = (LCM (Ci,...;Cm-1)). 

Proof of item 4. Since K is a submodule of S, so the generators 
Ki,...,Km-1 of K can be expressed as K; = Die Bp; for some Bx € R. 
Thus, K = ST where T = [Gji)i,j=1,..,m—1 © Mm-1,m-1: 

By Theorem 2.1, T = UDV for some invertible U,V € Mm-1m—1 and 


t 0 0 
0 te a 

a where ¢; are the elementary divisors of T. 
0 0 --: bm—1 | 


K=ST=SUDV <> Kv "=(SU)D. 


Since the columns of KV-! = [Ki,...,K’,_,] is another basis for K; and 
similarly, columns of SU = [p},...,p},-1] is another basis for S. Thus, we 
have the equality of the matrices 


Ga heal a hV = (SU) D= Tie teil: 
Thus consider the following map 


é: RxRx-:xR + (SU)/(KV™")3S8/K 


(T1,72,. . mi) > TP} dew an eb Tmt 
We see that 
m—1 
ker(¢) = {("1, wae tea) Ee Rm! | >- rip c Kv-} 
i=1 
= (t1) x (tg) mores Gai) 
Therefore, 


S/K © (SU) /(KV72) © R/ (ty) x B/(ta) x +++ x B/(tm—1)- 
C 


We want to emphasize that Theorem 3.1 is true under the condition 
that the polynomial ring is univariate, and not true in general for multi- 
variate polynomial rings. 

Next, we will continue with the notations used in Theorem 3.1, and 
show that SNF stated in Theorem 2.1 can be used to identify a some spe- 
cial properties of syzygy modules. 
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Theorem 3.2. Let S = [p,...,Dm—1] € Mmm-1 whose columns are the 
minimal set of generators of the syzygy module S of fi,...,fmj and K = 
[Ky,...,Km-1] whose columns are the minimal set of generators for the 
“Koszul-like” syzygy submodule K over the univariable polynomial ring R = 
K{2] obtained in Theorem 3.1. Let K =UDV be the SNF for K. Then the 
last row Um of the matriz U~* is a constant multiples of fi,...,fm- 


Proof. First, we express part 2 of Theorem 3.1 in term the following matrix 
equation: 


Cy Oy. Hes 0 
0 Cy: 0 
[Cipts+++1Cm—1Pm~-1) = [P1,-+- 5 Pm-1 os ; 
0 0 Cm-1 
= KQ for some Q € Mn-1m-1 
id 0 6 OU 
0 dy. te3-~ 00 
= UDVQ=U Ee | EQ 
0 0 > dm-1 
0 0 0 


where U,V aie invertible U € Maw V € Mm-i1m-1, D € Mmm-i- 


Multiplying both sides by U~! yields 


d, 0 0 

0 dy -: 0 
U-"Cyp1, Cope, ..+)CmtPm-1] = |i i} > | VQ, 

0 O 0 


and the last row of this matrix equation gives 
UZ (C1p1, Cope, seay Cin—1Pm-1| = (0, 0, cae , 0 


=> Un(Cipi) =Ci(Ump;) =0 => Ump,=0, Vi=1,...,m—-1 
=> Ur, is orthogonal to p;, Vi=1,...,m—1. 
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Since {p1, P2,.-+;Pm—1} is a minimal set of generators for the syzygy module 
of fi,..., fm, we must have [f,,..., fm] is orthogonal to each p; for + = 
1,...,m—1. Since S = [p,,..., Pm—1] and rank(S’) = m—1, the rank of null 
space of S is one. Hence U,, = 6[f1,..., fm| for some B € R = K{z]. Now, 
the fact that U is invertible over R implies that det(U) € K. Therefore, 
6 €K, otherwise det(U) would have a non-constant polynomial as a factor, 
contradicting to the fact that U is invertible. 
Thus, we conclude that the last row U,, of the matrix U~ is a constant 
multiples of fi,..., fm- 
C] 


We shall use the following example to illustrate the the results in The- 
orem 3.1 and Theorem 3.2. 


Example 3.3. Consider [fi, fo, f3] = [1,27 — 1,22 +1]. One can compute 
a minimal set of generators for the syzygy module Syz2(f1, fo, f3) is formed 
g?—1 2°41 
by the columns of the matrix S = [pi,p2] = | —1 0 
0 —l 
The submodule K is generated by the syzygies in the columns of the 


matrix 
w?—-1 2341 0 
(Ki, Ko, Ks] = —1 0 g?—a+1 . 
0 -1 -(#%-1) 


Since Kz = —(2? ~2 +1)K, 4+ (x—1)Ka, the submodule K is minimally 
generated by Ki, Ky. We will write K = (Ki, Ko]. We note that Ky = py, - 
and —27K,+ Ko = po, that is 


e—-1l £341 1 ee—-1 «+1 
—£ 
—1 ) f 1 | = —1 xr 
0 —1 0 —1 
[P1, Po]: 


Following the notation of Theorem 3.1, Ci = Cp =1, and therefore, 


I 


Ktal|y 


(K +S) = (LCM(C,,@:))=(1). => S2K. 
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, : 1-2] fi « 
This can also be verified as the following. Since 01 =o 4) 


g—-1 23+1 e—1 xt+1 = 
ae tet ie sea ue pee 
Q 1 


0 =] 0 —1 


[Ky, Ko] 


(p1; DoJ = ST where T = k i ‘ 


In terms of Theorem 3.1 (4.), the elementary divisors of T are: ty = tg = 1. 
Hence S/K = R/(1) x R/(1) = 0, that is, SSK. 
Now, compute the SNF of K 


z—1 23+1 1] fi Oo 10 
kK = UDV=|-1 0 01 |0 1 4 i where 
0 —l1 0; 10 0 


0 -!l 0 
Us S10. 10 -1 |. 
1 2-1 241 
The last row of the matriz U~+ is U3 = [1,27-1,2°+1] = [fi, fo, fal, which 
verifies the result of Theorem 3.2. 
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Abstract 


The central concept in the harmonic analysis of a compact group is the 
completeness of Peter-Weyl orthonormal basis as constructed from the matrix 
coefficients of a maximal set of irreducible unitary representations of the group, 
leading ultimately to the direct sum decomposition of its L?- space. A Peter- 
Weyl theory for a semicomplete orthonormal set is also pos- sible and is here 
developed in this paper for compact groups. Existence of semicomplete 
orthonormal sets on a compact group is proved by an explicit construction of 
the standard Riemann- Lebesgue semicomplete orthonormal set on the Torus, 
T. This approach gives an insight into the role played by the L?— space of a 
compact group, which is discovered to be just an example (indeed the largest 
example for every semicomplete orthonormal set) of what is called a prime- 
Parseval subspace, which we proved to be dense in the usual L?~ space, serves 
as the natural domain of the Fourier transform and breaks up into a direct-sum 
decomposition. This paper essentially gives the harmonic analysis of the prime- 
Parseval subspace of a compact group corresponding to any semicomplete 
orthonormal set, with an introduction to what is expected for all connected 
semisimple Lie groups through the notion of a K-semicomplete orthonormal 
set. 
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§1. Introduction. 

Harmonic analysis on a compact group is mainly a direct consequence of 
the famous Peter- Weyl theory which gives a consistent method, via the com- 
putation of the matrix coefficients of its irreducible unitary representations, of 
deriving a complete orthonormal set which is immediately responsible for the 
direct-sum decomposition of its L?— space and regular representation. Even 
though such a complete orthonormal set is non-existence for non-compact 
topological groups and hence the harmonic analysis on non-compact topolog- 
ical groups, as we know for connected nilpotent and semisimple Lie groups, 
has had to be developed through other means notably via the differential 
equations satisfied by the (spherical) functions derived as matrix coefficients 
of irreducible unitary representations constructed from parabolic and cohomo- 
logical inductions and the completeness afforded by the Plancherel theorem 
(which in the final analysis still depends on the availability and properties of 
the discrete series (known to be the irreducible unitary representations corre- 
sponding to some complete orthonormal set) of some distinguished compact 
subgroups), it still found to be appropriate (and to have a sense of finality) 
to have some forms of Peter-Weyl results on such non-compact topological 
groups. 

It is however possible to get at the decomposition of the regular represen- 
tation of a compact group G (for a start) via the indirect use of the notion of 
a semicomplete orthonormal set on such a group, leading to the consideration 
of a distinguished subspace of L*(G) which is established to be topologically 
dense. The study in this paper opens up this field of research by a detailed 
look at the compact case. The paper is arranged as follows. 

§2. contains a quick review of the well-known notion of a complete or- 
thonormal set on a compact group, giving the detailed of the aforementioned 
consistent way of constructing such a set through Peter- Weyl theorem which 
then leads to the direct-sum decomposition of its [?—space. The concept of a 
semicomplete orthonormal set on a compact group G is introduced in §3. with 
constructible examples (prominent among which is the Riemann-Lebesgue or- 
thonormal set)on the Torus, T, where we derived and used the properties of 
the Fourier and prime-Parseval subspaces of L?(G). Chief among these prop- 
erties is the topological denseness of every prime-Parseval subspace in L?(G). 
This takes us to the Fourter transform of the prime-Parseval subspace and its 
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direct-sum decomposition into invariant subspaces. The last section gives an 
introductory extension of the results of §3. on compact groups to connected 
semisimple Lie groups with finite center. 


§2. Fourier and Parseval subspaces for complete orthonormal set. 

A mutually orthonormal family {ya}ae4 in a Hilbert space, (H, (:,+)) is 
said to be complete (in H) if « € H is such that (z, ya) = 0 (for every a € A) 
implies z = 0. This means that a family {ya}aea of mutually orthonormal 
members of H is complete whenever it can be shown that the zero element 
of H is the only non-member of the family that is mutually orthonormal to 
all members of the said family. Two other equivalent methods of confirming 
the completeness of the family {xa}aea are as follows. 

2.1 Lemma. ([5.], p.3) Let {xahaea denote a mutually orthonormal 
family in a Hilbert space (H, (-,-)). The following are equivalent: 
(a) Every x € H can be expressed as © = Yyye4(X,Xa)Xa- 
(b) Buery « € H satisfies || x |P= ven | (es Xa) [- 
(c) {xa}aea ts complete in H.O 

The informed reader would observe that (a) of (2.1) is a Fourier series 
expansion of x while (b) of (2.1) is its Parseval equality, both with respect 
to {xa}aea. The import of this equivalence (in the light of (a) of (2.1) (re- 
spectively, (b) of (2.1))) is that every xz € H has a Fourier series expansion in 
terms of any known complete orthonormal family in H. We could then say 
that the subset H(x.) of H given as 


{ce Hers Sot, Xa)Xa; for some orthonormal family {xe}aea in H} 
cA 


(equivalently, the subset Hy(xq) of H given also as 


{x€H:|| a |?= Ss | (tz, Xa) |*, for some orthonormal family {y.}ae4 in H}) 
acA 


is exactly H if, and only if, {xa}aea is complete. Indeed another version 
of the equivalence of Lemma 2.1, whose formulation serves as our point of 
departure, is given as follows. 

2.2 Lemma. Let {xahaea denote a mutually orthonormal family in a 
Hilbert space (H, (-,-)). The following are equivalent: 
(a) H(Xa) =H 
(b) Hop(Xa) = H 
(c) {xa}aea ts complete in H.O0 
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2.3 Remarks. It may be safely conjectured that the Fourier subspace 
H(xXq) as well as the Parseval subspace Hy(xXo) (of a Hilbert space H) with 
respect to a complete mutually orthonormal family will always be equal to 
H, It will be a delight to study the disparity between the Fourier subspace 
H (Xa) as well as the Parseval subspace Hy(Xa) (of H with respect to the 
mutually orthonormal family {yo}ac4) and their inclusions in H, when the 
family {va}aea is not complete. 

For example, if the family {ya}ae4 of mutually orthonormal members 
in H is such that (x, x.) = 0 (for every a € A) does not necessarily imply 
whether x = 0 or x #0, it possible to then have that 


0 <|| 2 P= >_> | xa) P= 0, 


acA 


showing in this case (for the family {xo}oea in which (2, x.) = 0 (for every 
a € A) does not necessarily imply whether x = 0 or z # 0) that we now 
have Hu(xo) = {0} (= H(xa) # H, showing that both subspaces are too 
small and far from being equal to H). This shows at a glance the importance 
of completeness of the family {va}ae4 in the consideration of the Parseval 
equality, for the non-triviality of these two subspaces H(xq) and Hy(xq) and 
for the sustenance of the relationship of equality (of Lemma 2.2) between 
(Xa) and Hy(xo). 0 

However, and as it shall be shown in the next section, these two sub- 
spaces, H(x.) and Hy(xa) may be considered for an appropriately chosen 
not-necessarily complete orthonormal family {xq}ae4 and with which they 
would still be found not to be too small in sizes (in comparison with H). 
This choice of a not-necessarily complete orthonormal family {ya}ae4 would 
equally help and be appropriate in order that both H(xvq) and Hy(xa) be 
lifted to all of H. All this in a moment. 

A well-known method of computing complete orthonormal family of func- 
tions is via the matrix coefficients of irreducible unitary representations of 
a compact groups G which is then used to decompose L?(G) into invariant 
subspaces, leading to the decomposition of the right regular representation 
on G (which sadly, does not generalize to non-compact topological groups). 
Here is the technique. a 

Denote the dual of a compact group G by G, consisting of all its equiv- 
alence classes of irreducible unitary representations. For A € G denote by 
ui the corresponding matrix coefficient representative of the class 1 whose 


Bess am 
degree is also denoted by d(\). Then the set 


{dau : +€G,1<i,j < dQ} 


consists of a maximal set of complete orthonormal family of functions in 
L*(G) and (hence) every f € L?(G) can be expanded as 


d(d) 


f= dod ) Dh os) uy 


EG 
(with convergence in the norm of Aes (G)) whose Fourier transform 
f:@ > May(C): 04 fA) = (FWD, 
is given as FO)y ae ux) (where Ma)(C) denotes the algebra of matrices 


with entries in C and degree d(X)). It then follows that for any compact 
group G, the Fourier subspace L*(G)(1/d(A)uj;) of L?(G) is given as 


d(\) 
2(a\(/dA)ud) = {fe PG): f= > od) ,uj)uj} = PG) 
AEG aj 


‘(=L*(G)p(/d(A)ujs), the Parseval subspace of L*(G)), with respect to the 
family {\/d(\)us: d € G, 1 < i,j < d(A)}. We then have the abstract 
direct-sum decomposition of L?(G) given as 


aA) 


-OO" 


where HA := ye Cu. This is the content of Peter-Weyl Theorem, [5.], 
and we shall refer to the set 


{JdA)jur: +164, 1<4,7 <d(d)} 


as the standard Peter-Weyl orthonormal set on G. 

The inability of being able to get an orthonormal family in L*(G) for a 
non-compact topological group G in the above tradition of Peter-Wey] is the 
first stumbling block to harmonic analysis on such groups, which has been 
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considerably understeod and completely developed via a rigorous treatment 
of the rich structure of differential equations satisfied by matrix-coefficients 
of members of each of the classes in G, [2]. This paper presents a constructive 
method of getting a not-necessarily complete orthonormal set which is close 
enough to being a complete orthonormal family in an arbitrary Hilbert space 
(H, (:,-)) and/or in L(G), for a compact group (and introduced the same 
technique for a semisimple Lie group) offering a more general Fourier series 
expansion of each member of an appropriate subspace of H and/or L(G). 

Starting with a compact group (before extending the notion to all con- 
nected semisimple Lie groups, with finite center, via its Iwasawa decomposi- 
tion) we would however not approach harmonic analysis on the groups via 
the completeness (and consequent denseness) of the standard Peter-Weyl or- 
thonormal set, but via a denseness in the L?—space which would be found 
to be possible from an almost complete orthonormal set. 


§3. Semicomplete orthonormal set in a compact group. 

The existence of different special functions and polynomials of math- 
ematical physics, which have been established to be orthonormal in vari- 
ous semisimple Lie groups (compact and non-compact types), is well-known. 
However the absence of completeness of these orthornormal families (under 
the structure of their individual corresponding groups) is the first stumbling 
block to a direct Peter-Weyl harmonic analysis of them. In this section we 
shall define and study the concept of a semicomplete orthonormal family in 
a compact group in order to extend this concept to the harmonic analysis of 
all semisimple Lie groups in the next section. 

3.1 Definition. (Semicomplete orthonormal family) Let G denote a com- 
pact group and let the members of the non-empty set A be ordered such that 
A = {aj}ij. An orthonormal family {xi hedca L?(G) is said to be semi- 
complete if given € > 0 there exist some non-zero scalars 


Viot*' Veo »Pry**' By EC 


andn €N such that 


d(A) n n 
| S>d(d) Soy, Uy )un — Soy S> Bis ys (F, Kod Ned lla< ¢ 
1G tg j=l i=l dea 


for every f € L°(G).0 


90: 


The quantity 


d(A) 
ddO) DCP uaa 
AEG ijel 


in Definition 3.1 above may be replaced with f (due to the Peter-Weyl The- 
orem), so that the other quantity 


ay By Do (Fs Kat Xes 


t=1 ob EA 


(in the same Definition above) should be seen as the total contribution of 
{Ket bade , in L*(G) in its bid to attain f, Thus the informed reader would 
see that the inequality in na * above simply gives a measure of how 
close to the completeness (of {/d( d)ui}) ) is the orthonormal set {x as od eA: 


The standard Peter-Weyl orthonormal basis {./d(A)u; iy} used in the above 
Definition 3.1 may be replaced by any other known complet orthonormal 
set {vy},ep in L?(G) while the concept of a semicomplete orthonormal set 
(for {xa} ica) could also be defined for an arbitrary Hilbert space, H, so 
as to have what may be generally called a semicomplete orthonormal set in 
H with respect to (the complete orthonormal set) {vy},epn in H. If in this 
general case the set {v,},ce in H is also not necessarily complete, we may 
arrive at the notion of a relative semicomplete orthonormal set for {x ,;} aleA 
in H with respect to {v,},ee8 in H. Thus Definition 3.1 may therefore be 
seen as giving semicompleteness of Ned tal jeq im L?(G) with respect to the 


standard Peter-Weyl orthonormal basis { d(r)uss }. 

It is clear that every complete orthonormal set in L?(G) (or i in any Hilbert 
space, /7) is automatically semicomplete; simply choose A = G 1; = By = 1, 
but not conversely. An inductive method of immediately constructing a 
semicomplete orthonormal set in a compact group is by a method of selective 
omission of some number of members in any known complete (or of the 
standard Peter-Weyl) orthonormal set with a controlled bound. The control 
of the bound in the method of selective omission would be achieved using 
the Riemann-Lebesgue Lemma. 

This method, as contained in the following, equally gives an existence 
argument for the concept of a semicomplete orthonormal set in a compact 
group. 


JBie 


3.2 Lemma. (Existence of a semicomplete orthonormal set: the standard 


Riemann-Lebesgue orthonormal set on the Torus, T) There exist ro € T for 
which F 
N 
| (fier) I< d(Xo)?’ 


for every f € L?(T), |X |>| Ao | and 1 < k,m < d(Ao). Moreover, 


{V/d(\ur : 4 € G\ fro}, 14,9 < d(A)} 


is a semicomplete orthonormal set on T. 
Proof. Since the dual group T is discrete, so that 


lim (f, oy) = lim FOX Jig =O (by the Riemann-Lebesgue Lemma), 


|A\-too |Al00 


it follows that there are (infinitely) many possible  € é (choose such one 
do) with | A [>| Ao | for which | (f, um) =| (f, Uda) — 0 I< Woy for every 
f € L°(T) and 1<k,m < d(Ao), as required. Hence, 


d(\) d(A) Aro) 


I Soa) SO (feud yud— SS dQ) SO (f,ud)ud llo=ll do) SO (f, uleyady 


AeF i=l AEF\{ro} =I wet 


< d(Ao) ped, | (f, ui?) |< e, for every f € L*(T).0 

The technique ah ‘Lemma 3.2 may be extended as ee oe 
choose (as assured by the Riemann-Lebesgue Lemma) MY on ,: € T for 
which 


ee ee 
do | fii) I< swe 


where | A |> max{| Ms” |, NM) |, .--} and f € L*(T). Then, with proof 
essentially the same as in Lemma 3.2, the set 


{Jd\)ud: € G\ OY AP}, 147 < day} 


is a semicomplete orthonormal set on T. We shall henceforth refer to the 
semicomplete orthonormal set 


{JdA)jud : 6 B\ OP AY, 145 < day} 


lla 
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as the standard Riemann-Lebesgue (semicomplete) orthonormal set on T (be- 
ing in correspondence with the standard Peter-Weyl (complete) orthonormal 
set, {x/d(A)uj}.) 

Other non-standard examples of Definition 3.1 may be deduced from the 
numerous special functions of mathematical physics where their correspond- 
ing non-zero scalars y; and @,; in Definition 3.1 could be calculated from. 

3.3 Remarks. In contrast to the zero-subspace Hy(xq) of Remarks 2.3 
we may, in the context of a semicomplete orthonormal set {xa}eea in a 
Hilbert space (H, (-,-)), consider the subspace 


Ay (Xe) = {t € H: (€,xXa) =0, (for every a € A) implies x = 0}, 


for some orthonormal set {ya}aca in H. It is clear (from Lemma 2.2) that 
Hyy(Xa) = H (hence equal to H(xa) and Hy(Xa)) if, and only if, {va}aea is 
complete in H and that, when {xa}aea is semicomplete in H or in L7(G), 
both Ay(xa) and Hy(Xq) are non-zero: an example may be seen from using 
the standard Riemann-Lebesgue orthonormal set on T. In general, we have 
the following. 

3.4 Lemma. Let (H, (-,-)) denote any Hilbert space. Then 


(Xe) G Fy (Xe) 


for any semicomplete orthonormal set {yasaca in H. 
Proof. Choose any « € H(yq), then & = )oyea(®)Xa)Xa. Now if 
(XL, Xa) = 0, for every a € A, then 


c= Sa, Xe) Xa = S°(0)xXe = 0; 


acA acA 


showing that z = 0 as required. UO 

We shall refer to Hy,(Xa) as the prime-Parseval subspace of H and the 
choice of this term is further reinforced by the following facts. 

3.5 Lemma. (cf, Lemma 2.2) Let {xa}aca denote a semicomplete or- 
thonormal set in a Hilbert space (H,(-,-)) and let x € H. Then x € Hy(Xq) 
whenever || x |P= dooce, | (2) Xe) [? - | 

Proof. If || x ||?= vaca | (©) Xa) |? and | (2, Xa) |= 0 (for every a € A), 
then || x |?= ooc4 | (ts Xe) P= Yeea(0) = 0; showing that « = 0. Hence 
x € Hy (Xa). 

Lemma 3.5 shows the first partial connection between the satisfaction 
of Parseval equality, on one hand, and membership in the prime-Parseval 
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subspace, on the other. The last Lemma may also be seen as saying that the 
subset of H given as 


{cs € HH: || a |P= > | (x, Xo) |?, for any orthonormal set {Xahoca} 
acA 


is also a subset of Hiy(xq), with clear equality when {xa}aeA is complete. It 
will be satisfying to also have the reverse inclusion, 


Hy (Xe)  {e € H: || x |P= > | (w, xa) |?, for any orthonormal set {xo }aca} 
acA 


due to the importance of the Parseval equality in the fine properties of Fourier 
transform. We shall deal with this concern in Lemma 3.12. 

Even though a semicomplete orthonormal set {Xa}qea in L*(G) (or ina 
Hilbert space (H, (,-, ))) may not be dense, as it is generally expected of a 
complete orthonormal set, we may still however employ this orthonormal set 
to construct some dense subspaces of L*(G) (or of a Hilbert space (H, (,-, ))) 
as follows. Indeed, the following results on the Fourier subspace for L?(G) 
are also valid for an arbitrary Hilbert space, (H,(,:,)) and for a relative 
semicomplete orthonormal set in H. 

3.6 Theorem. Let G denote a compact and let Wed tade 4 denote a 


semicomplete orthonormal set on G. Then L?(G) (Xai) is topologically dense 
in L?(G). 
Proof. Since every f € L?(G) may be expanded as 


d(A) 
S- dr) So (f,ud)us, 
1G tl 


(with convergence in the norm of L?(G)) it follows that for « > 0 we have 


a(A) 
E 
lf — Dla) DoF miiyuds lle< 5; 
AEG igj=1 


Now 


d(A) 


S29 2B YF ef) Xe US I — OQ) DOA) Il 
j=1 


j= i=l aJEA AEG i,j=1 
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d(X) 
+ |] SodQ) Sof, wus - a Yh = (f; Xed)Xed lla< 5+ 5 =«0 
eG ij=l aleA 


In more specific terms we have the following. 
3.7 Corollary. Let G denote a compact group and let {x,,} ica denote a 


semicomplete orthonormal set on G. Then every f € L?(G) can be expanded 


as 
f= uy a ye. (f; Xod)Xod 
j=l i=1 od cA 

for some 1;, Bij € C with convergence in the norm on L*(G). 0 

We may refer to the expansion of f in Corollary 3.7 as a semi-Fourier 
series expansion for f € L?(G) or H with respect to Xba a- A stronger 
form of Theorem 3.6 carved in the form of the equivalence of Lemma 2.2 
and which generalizes the fact that a mutually orthonormal family {xa}ae 
is complete (in a Hilbert space (H, (-,-))) if, and only if, H(xo) = H (ef. 
Lemma 2.2) is also possible when the mutually orthonormal family {va}aca 
is semicomplete in H. We prove this below in the special case of H = L(G). 

3.8 Theorem. Let G denote a compact group and let {Xai} oleA de- 
note a mutually orthonormal set on G whose Fourier subspace is denoted as 
L*(G)(Xqj)- Then L*(G)(Xqj) is topologically dense in L?(G) if, and only if, 
(ee, 4 18 semicomplete. 

Proof. That L? (G)(Xeg') is topologically dense in L7(G) if {Xai dedea is 
semicomplete is the content of Theorem 3.6. Now choose f € L? (G), then 

d(A) 


I Day 52( (f, udyur SS hy dof Xe8) xed lle 
AEG ij=l j=l i=1 od CA 
d(A) n n 
SI D0 4) SOF whut —F lle + WP — Dow D0 Ba SOA Xe) Xe Ile 
NEG 4j=l j=l i=1 aJeA 


ae a e(using the Peter-Weyl theorem and Corollary 3.7, respectively). 0 
This Theorem would enable us to see the Peter-Weyl series expansion of 
every f € L7(G), given as 
d(A) 


f= > odd) do (f,uh)ug 


‘eG tel 
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(with convergence in the L?—norm), as the restriction of the semi-Fourier 


series eEZTpANsion 
f= Lado Do (f; X0g)X 


j=l i=l ay JEA 


to the standard Peter- Weyl (complete) mutually orthonormal set {./d(A)u;s}. 
Indeed Theorem 3.8 leads to the same conclusion for the prime-Parseval sub- 
space L?(G)p (X04): 

3.9 Corollary. Let G denote a compact group and let {Xas} ated denote 
a mutually orthonormal set on G. Then L?(G)ry (Xej) 48 topologically dense 
in L?(G) if, and only if, {Xai haiea 48 semicomplete. 

Proof. Consider Lemma 3.4 in the light of Theorem 3.8. 0 

The inclusion L?(G)(x,3) € L?(G)i (x4, i) of Lemma, 3.4, when combined 
with both Theorem 3.7 and Corollary 3.9, implies the following. 

3.10 Corollary. L?(G)(x Xej) is topologically dense in L?(G)is(X,9)- 0 

The converse of Lemma 3.5 is now immediate for both L?(G)qy (Kad ;) and 
(even) Hj,(X,3) in any arbitrary Hilbert space, (H, (-,:)). 

3.11 Lemma. (cf Lemma 2.2) Let G denote a compact group and let 
{Xa}aca denote a mutually orthonormal set on G. Then f € L?(G)ig(Xa) tf, 
and only #f, |l f llB= doaea | (f:Xa) |? - 

Proof. Let f € L?(G)x(xa). We may take f € L?(G)(xXa) due to Corol- 
lary 3.10; so that f = eet (f;Xa)Xa- Hence 


O=|| f — SOF, xa) Xe =H £2 — DO | (A Xe), 
aca aca 
as required. O 
Hence, the prime-Parseval subspace L?(G)y(x od) may finally be seen (for 
some orthonormal set Ugh ale A) 28 


TP(G)g(xq4) = Lf € L7G): WF B= D2 | (Fs Xag) P 


ajcA 


We now have enough preparation to introduce a Fourier transform f '-> f 
on the prime-Parseval subspace, LP? (G)ip (Xai): 


~ 


Consider f € L*(G) and for every a € A define the matrix f(a) whose 
entries are given as 


~ 


Fla)iz = fled). 


2352 


That is, Fladij = (f,x,4), for 1 < i,7 < n. The Parseval inequality of 
L?(G)5p(X a3) (in Lemma 3.11) therefore becomes || f |]2= Scaca || f(a) ll?, 
for every f € L?(G)i(x,3), where || f(a) |? is the Hilbert-Schimdt norm of 


the matrix - 
Fla) = (f@)ahjar = Fe) ja 


In other words, and in terms of our choice of indexing A, we have 


Il f lb= ys >= I Fed) IP 


i,jal od cA 


for f L(G) (x44) 

3.12 Definition. Set L(A) as the space of matrix-valued functions yp 
on A with values in UP2_, Mn(C) satisfying 
(i) p(od) € M,(C al € A and 


(it) Sea 1 aleA | y(a wales oo. 
The inner product ie ; on L?(A) given as 


(y,)) = 3 S> ir(y )*), 


ij=1 oJeA 


vy, € L*(A) converts (£?(A),(-,-)) into a Hilbert space. We can then estab- 
lish a connection between the prime-Parseval subspace L?(G)ry ea) (which 
is a Hilbert subspace of L*(G)) and L*(A). 

3.13 Theorem. (Fourier image of the prime-Parseval subspace) Let G 
denote a compact group and let {Xoi bade 4 denote a semicomplete mutually 
orthonormal set on G. Then the map j 


H.: L?(G\ip(xq3) + 17(A): f HUGS) = F 


is an isometry of L?(G)ig (Xd) onto L?(A).0 

Theorem 3.13 is very familiar when the semicomplete a ae orthonor- 
mal set {x aj} ajca is the complete mutually orthonormal set {/d(A)us}. We 
do not yet know the general connection between the set A and the dual group 
G, except in the special cases of the standard Riemann-Lebesgue (semicom- 
plete) orthonormal sets on T. We however see A as a general form of G which 
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may take the usual form of Gin specific cases. If we set 


Hf = S0Cxe4; 


j=l 


fora = a} € Aandi € {1,-+- ,n}, then the Hilbert subspace L? (Gin (Xed) 
of L?(G) has the direct-sum decomposition 


L*(G)ip(xa) = DD HF. 


acA i=l 


The results of this section laid a foundation for harmonic analysis of the 
prime-Parseval subspace Fg (Xoi) with respect to a semicomplete orthonor- 
mal set {Xai } ale A in a Hilbert space, H. Having considered the case of the 
Hilbert space L?(G), for a compact group G, in this section it will a delight 
to use these foundational results (on both Ai,(x ad) and L?(G)in(x 3) in the 
understanding of further properties of D?(G)in(X 9) in the full sight of the 
semicompleteness of {x e) sie: We shall give a very short introduction to 
this type of study for a connected semisimple Lie group in the next section. 

It is clear from Lemma 3.2, for standard (Riemann-Lebesgue) examples of 
a semicomplete orthonormal set in an arbitrary Hilbert space (H, (:,-)) or in 
L*(G), that the non-zero constants y; and /;; would always be y; = Bj; = 1 
for 1 <i,7 <| GX pu re) --»} |. However, for non-standard examples of 
a semicomplete orthonormal set in an arbitrary Hilbert space (H, (-,-)) or in 
L?(G), the semi-Fourier series expansion of Corollary 3,7 may have to be 
broken down in order for general expressions for 7; and {j; to be known. A 
first result along this line is the following. 

3.14 Lemma. Let {Xoi hate 4 denote a semicomplete orthonormal set in 
a Hilbert space (H,(-,-)) and let x € H. Then 


(©, Xai) = KP(®, Xai)» 


for1<i<n. In particular, ¥;6% = 1. 

Proof. We have that (x, Xat,) = ae i ei Bij DUeA (x, Xed) (Xed> Xat,) . 
Due to the orthogonality of the set {x,;},7¢4 the above equality reduces to 
(2, Xai) = ViPis(C Xos), for l Si <n as required. 

Now (1 — ¥i4i) (x, Xqi) = 0 from where we have ¥;4; = 1.0 
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§4. K-semicomplete orthonormal set in a semisimple Lie group. 

The success in §3. of the use of the notion of a semicomplete orthonormal 
set in the harmonic analysis of a compact group, culminating in the extraction 
and elucidation of the prime-Parseval subspace as well as its Fourier image, 
shows the central importance and the correct us of Parseval equality and 
the concept of completeness (of an orthonormal set) in the abstract Peter- 
Weyl] theory of a compact group and in the understanding of the hitherto 
unknown subspaces of L?(G) under the influence of the Fourier transform. 
This study (which led us to the consideration of the prime-Parseval subspace 
L*(G ip (Xe) corresponding to a semicomplete orthonormal set {v of Nag ea on 
G) is reminiscence of and may be compared with the extraction and harmonic 
analysis of the Schwartz algebra in the L?—theory of semisimple Lie groups 
which was started in the Yale thesis [1(a.)] of James Arthur (continued and 
completed in two later manuscripts, [1(b.)] and [1(e.)]). In a more recent 
publication, harmonic analysis of other spaces of functions on semisimple Lie 
groups, namely of the space of spherical convolutions, has been introduced 
in [3.] leading to the explicit construction of the corresponding Plancherel 
formula for such functions. The present paper has also introduced the Fourier 
and prime-Parseval subspaces of L(G) (or of any arbitrary Hilbert space, 

Having shown in §3. the essential importance of the Parseval equality 
(which is the precursor of the Plancherel formula) in the consideration of the 
actual subspace of L?(G) under the natural action of the Fourier transform, 
we shall here consider studying the same theory (of a semicomplete orthonor- 
mal set) but for all semisimple Lie groups, having removed the impediments 
posed by the completeness for orthonormal sets on such Lie groups. 

It is well-known that orthonormal sets (of functions and polynomials) are 
numerous and readily available in the L?—space (and more recently in some 
distinguished subspaces of the L?"—spaces [4.]) of semisimple Lie groups. 
Indeed every semisimple Lie group has its corresponding orthonormal set, an 
example is G = SL(2,R) and its Legendre functions. 

Even though these sets of orthonormal functions and polynomials are 
central to harmonic analysis on these groups, their direct importance in or 
contribution to the decomposition of (sub-)spaces of L*(G) or expansion of 
their members is not yet-known. In the outlook of the present section (and 
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of the entire paper) any orthonormal set on a semisimple Lie group known 
to have been K—semicomplete (in the sense to be soon made precise) could 
be a basis of some subspaces of L?(G). 

4.1 Definition. (K—semicomplete orthonormal set) Let G= KAN de- 
note the Iwasawa decomposition of a connected semisimple Lie group G with 
finite center. An orthonormal set {Xa}aca onG is said to be K—semicomplete 
whenever its restriction to K, written as {(Xe)|c tac 18 @ semicomplete or- 
thonormal set in L?(K).0 

It is relatively easy to construct a K—semicomplete orthonormal set on 
any connected semisimple Lie group G, from any given semicomplete or- 
thonormal set on K as follows. 

4.2 An example, Choose any of the numerous orthonormal sets {£,}aeA 
in L*(K) as constructed in §3. and, for every = kan € G, define the map 
Xai: G—oCas 

Xa(z) = Xa(kan) = eflonle (k), 
where f : AN — C satisfies 
(¢) f(1) =0, 
(i4) fay eA") dadn = 1 and 
(itz) San g (Kan) (efr)+F(e))dadn = g(k), forg € L?(G),a € Am EN 
and the normalized Haar measures da and dn on A and N, respectively. 
Proof. Observe that since 


Xa(z) = Xa(kan) = efor) €.,(k), 
then for any k € K 
Xa(k) a Xalk ot 1) = ef ONE, (k) = Ea(k). 


For any 01,2 € A, we have 


Oar» Xea) = im e2(I(0")) dad) a, (k)Eca(B)dk = (as Eon) 


and 
I xa l= f Cf eR" dad) | fas () P dk =x IB 2 
K JAN 
showing that {xa}ae, is an orthonormal set on G. Its K —semicompleteness 
is also shown as follows. For a pre-assigned € > 0, we have that 
d(A) 


I do 4A) So (a, ud ud — Soy 2 Big DO 9 X at) Xee lle 
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~ afd) a 
=I 40) Puhr - [ Ef alean) eM) dan}: 
AEG ij=l K JAN 


ar > Bij > €,i(k) dk bod lle 


j=l =1 dea 
d(A) 


=I] 54a) S (9, wdyurs — Soy Bs So Gs basa Ila O 


neG ij=l j=l i=l dea 


For any K—semicomplete orthonormal set {ya}aea on G the correspond- 
ing Fourier subspace L?(G)(x.) of L?(G) is also given as 


P(G)(xa) ={f EPG): f= Df Xa)xah 


acA 


while the prime-Parseval subspace is 
LD? (G)a (Xe) ={f € 17(G): (f,xa) =0 (for every a € A) implies f = 0}. 


Clearly L*(K)ig(/d(A)us,) = L?(K) (from Lemma 2.2 (#4), both sub- 
spaces L*(K)(Xq) and L*(K)i,(X) are topologically dense in L?(K) (from 
Theorems 3.6 and 3.8 and Corollary 3.9) and there exists an isometry of 
L?(K)ig(Xa) onto L?(A) (from Theorem 3.13). We shall resume the study 
of the subspaces L*(G)(Xq) and L?(G)y(Xq) (for connected semisimple Lie 
groups, G) in another paper. 
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Abstract 


We give the exact contributions of Harish-Chandra transform, 
(Hf)(A), of Schwartz functions f to the harmonic analysis of 
spherical convolutions and the corresponding L?— Schwartz al- 
gebras on a connected semisimple Lie group G (with finite cen- 
ter). One of our major results gives the proof of how the Trombi- 
Varadarajan Theorem enters into the spherical convolution trans- 
form of L?— Schwartz functions and the generalization of this 
Theorem under the full spherical convolution map. 
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1 Introduction 


Let G be a connected semisimple Lie group with finite center, and denote 
the Harish-Chandra-type Schwartz spaces of functions on G by C?(G), 0 < 
p <2. We know that C?(G) Cc I?(G) for every such p, and if K is a maximal 
compact subgroup of G such that C?(G//K) represents the subspace of C?(G) 
consisting of the K—bi-invariant functions, Trombi and Varadarajan ((9.]) 
have shown that the spherical Fourier transform f + fi is a linear topological 
isomorphism of C?(G//K) onto the spaces Z(%°), € = (2/p) — 1, consisting 
of rapidly decreasing functions on certain sets §* of elemetitary spherical 
functions. 

We show the existence of a hyper-function on both G and $* (here named 
a spherical convolution) whose restriction to the group identity element, e, 
coincides with the spherical Fourier transforms, f > f, , of Schwartz functions 
f on G and which affords us the opportunity of embarking on a more inclu- 
sive harmonic analysis on G. Indeed [8.] contains a more general Plancherel 
formula for the collection of these functions. As a function on G its series 
expansion is in the present paper studied. We show that, aside from the 
fact that the spherical Fourier transforms, f FO), is the constant term of this 
series expansion, there is a region in G where the spherical convolution is 
essentially f(A). Various algebras of these functions are thus studied and 
ultimately embedded in £?(G). It is however clear that the results in [8.] and 
in the present paper may be extended to include what may be termed as the 
Harish-Chandra-type Schwartz spaces of Eisenstein Integrals on G. 


The following is the breakdown of each of the remaining sections of the 
paper. §2. contains the preliminaries to the research containing the struc- 
ture theory, spherical functions and Schwartz algebras on G, while the series 
analysis of spherical convolutions on G is the subject of §3, where we also 
extend the Trombi- Varadarajan Theorem to all spherical convolutions. The 
relationship existing among the Schwartz algebras of functions and those of 
spherical convolutions is considered in §4. 


2 Preliminaries 


For the connected semisimple Lie group G with finite center, we denote 
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its Lie algebra by g whose Cartan decomposition is given as g = t@p. Denote 
by @ the Cartan involution on g whose collection of fixed points is t. We also 
denote by K the analytic subgroup of G with Lie algebra t. K is then a 
maximal compact subgroup of G. Choose a maximal abelian subspace a of p 
with algebraic dual a* and set A = expa. For every \ € a* put 


gv. = {X €g: [A,X] = \(4)X,VEH € a}, 


and call \ a restricted root of (g,a) whenever g, 4 {0}. 


Denote by a’ the open subset of a where all restricted roots are # 0, and 
call its connected components the Weyl chambers. Let at be one of the Weyl 
chambers, define the restricted root positive whenever it is positive on at 
and denote by A?* the set of all restricted positive roots. Members of At 
which form a basis for A and can not be written as a linear combination of 
other members of A* are called simple. We then have the Iwasawa decom- 
position G = KAN, where N is the analytic subgroup of G corresponding 
to N = yeas ga, and the polar decomposition G = K - cl(A*t)-+ K, with 
A* = expat, and cl(At) denoting the closure of A*. 


If we set M = {k € K : Ad(k)H = H, H € a} and M' = {k € 
K : Ad(k)a C a} and call them the centralizer and normalizer of a in K, 
respectively, then (see [5.], p. 284); (i) M@ and M’ are compact and have 
the same Lie algebra and (ii) the factor = M’/M is a finite group called 
the Weyl group. to acts on ag as a group of linear transformations by the 
requirement 

(sd)(H) = Ns“), 
H Ea, s € tw, » € ag, the complexification of a*. We then have the Bruhat 
decomposition 
G=| | Bm.B 
s€to 

where B = MAN is a closed subgroup of G and m, € M' is the represen- 
tative of s (ie, s = m,M). The Weyl group invariant members of a space 
shall be denoted by the superscript ® while | tv | represents the cardinality 
of tv. 


Some of the most important functions on G are the spherical functions 
which we now discuss as follows, A non-zero continuous function y on G shall 
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be called a (zorial) spherical function whenever y(e) = 1, y € C(G//K) := 
{g € C(G): g(kivke) = g(x), ki, ko € K, x € G} and fey = (f *y)(e)+¢ for 
every f € C.(G//K), where (f*g)(x) = fo f(y)g(y7'z)dy. This leads to the 
existence of a homomorphism 4 : C,(G//K) -> C given as X(f) = (f *y)(e). 
This definition is equivalent to the satisfaction of the functional relation 


a plky)dk = o(a)yl), 2,y eC. 


It has been shown by Harish-Chandra [(6.] that spherical functions on G 
can be parametrized by members of ag. Indeed every spherical function on 
G is of the form 


~r(x) = i eP-PH (ek) dk Dd € at, 
K 


2 = § Dd yeat M+ A, where my = dim(ga), and that yp, = y, iff \ = su for 
some s € ty. Some of the well-known properties of spherical functions are 
p-rxlo) = gr(z), p-a(z) = Ox(2), | ale) |S perl), | pale) IS vier(a), 
p-ip(t) = 1, X € ag, while | ~)(x) |< Yo(z), A € ia*, c € G. Also if 2 is the 
Casimir operator on G then 


Oy, = —((A, A) + (P,P) Pars 


where A € ag and (A, u) := tr(adH adH,) for elements H), H, € a. This 
differential equation may be written simply as Ny, = 7(Q)(A)y~), where 
A ++ ¥(2)(A) is the well-known Harish-Chandra homomorphism. The ele- 
ments H), H, € a are uniquely defined by the requirement that A(H) = 
tr(adH adH,) and u(H) = tr(adH adH,) for every H € a (([5.], Theorem 
4,2). Clearly Qy9 = 0. : 


Due to a hint dropped by Dixmier [4.] (cf. [9.]) in his discussion of some 
functional calculus, it is necessary to recall the notion of a ‘positive-definite’ 
function and then discuss the situation for positive-definite spherical func- 
tions. We call a continuous function f : G > C (algebraically) positive- 
definite whenever, for all 71,...,2%m in G and all a1,...,@m in C, we have 


St and; f (ap x5) > 0. 


i,j=1 
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It can be shown (cf. [5.]) that f(e) > 0 and |f(z)| < f(e) for every cE G 


implying that the space P of all positive-definite spherical functions on G is 
a subset of the space $! of all bounded spherical functions on G. 


We know, by the Helgason-Johnson theorem ([7.]), that 
Fi =a°t+iC, 


where C, is the convex hull of {sp : s € tw} in a*. Defining the involution 
f* of f as f*(z) = f(x), it follows that f = f* for every f € P, and if 
( € P, then \ and A are Weyl! group conjugate, leading to a realization of 
P as a subset of tv \ a%. P becomes a locally compact Hausdorff space when 
endowed with the weak *—topology as a subset of L(G). 


Let 
Yo(z) = [| exp(-p(H(ok)))ak 


be denoted as =(z) and define 0: G+ C as 
o(x) = ||X|| 


for every x = kexpX €G, k € K, X € a, where || - || is a norm on the 
finite-dimensional space a. These two functions are spherical functions on G 
and there exist numbers c, d such that 


1 < E(a)e8% < c(1 + a(a))?. 


Also there exists r > 0 such that c =: f,8(x)?(1 + 0(x))"dx < oo (11, p. 
231). For each 0 < p < 2 define C?(G) to be the set consisting of functions 
f in C®(G) for which 


Ma,by (f) = supllf (a; 2; b)|E(x)-7/?(1 + o(2))"] < 00 


where a,b € MU(gc), the universal enveloping algebra of gc, r € Zt,x € G, 
f(a;b) := 4| 6 f(x - (exptb)) and f(a;z) = $|,_, f((expta) +x). We call 
C(G) the Schwartz space on G for each 0 < p < 2 and note that C?(G) is the 
well-known (see [1.]) Harish-Chandra space of rapidly decreasing functions on 
G. The inclusions 


CX(G) c C?(G) C IP(G) 
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hold and with dense images. It also follows that C?(G) C C%(G) whenever 
0<p<q< 2. Each C?(G) is closed under involution and the convolution, *. 
Indeed C?(G) is a Fréchet algebra ([10.], p. 69). We endow C?(G//K) with 
the relative topology as a subset of C?(G). 


We shall say a function f on G satisfies a general strong inequality if for 
any r > 0 there is a constant c = c, > 0 such that 


| fy) |< eq By 2) +o(y 2)" Va,yeG. 


We observe that if s = e then, using the fact that E(y-') = E(y) and 
o(y~') =a(y), Vy € G, such a function satisfies 


I fy) IS BY) + oy)" = GE) +o)", Vy EG, 


showing that a function on G which satisfies a general strong inequality 
satisfies in particular a strong inequality (in the classical sense of Harish- 
Chandra, [11.]). Members of C?(G) =: C(G) are those functions f on G for 
which f (91; °; 92) satisfies the strong inequality, for all g1, g2 € L (gc). We may 
then define C)(G) to be those functions f on G for which f(g1; +; go) satisfies 
the general strong inequality, for all 91, 92 € £ (gc)and a fixed x € G. It is 
clear that C (G)= C(G) and that U,¢qC™ (G),which contains C(G), may 
be given an inductive limit topology. The seminorms defining this topology 
will be explicitly given in §4. 


For any measurable function f on G we define the spherical Fourier trans- 
form f as 


7) = Ee f(x) p-a(w)dey, 


A € ag. It is known (see [3.]) that for f,g € L'(G) we have: 


(i.) (f xg) = f -G on ¥ whenever f (or g) is right - (or left-) K-invariant; 


(ii.) (f*)(~) = F(v*), 9 € &; hence (f*)* = Ff on P : and, if we define 
f#(9) = Speer f (kivke)dkidke,x € G, then 


(iié.) (f#)* = f on F. 
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We shall denote the spherical Fourier transform f(A) of f € C (G) by 
(Hf)(A) and refer to it as the Harish-Chandra transforms of f. Its major 
properties are well-known and may be found in [9.]. It should be noted that 
(A)(A) = fA) = Se f@)e-a)ay = Je Fw)ealy™ dy = Sg fealty e)dy 
= (f *y,)(e). That is, the Harish-Chandra transforms of f is the restriction 
of the function 

or+ (f * pr)(x) = 8,7 (x) 
on G to the identity element. It is therefore worthwhile to explore s),¢(z) 
in some details for all 2 € G in order to put its behaviour at z = e (as the 
Harish-Chandra transforms of f) in a proper and larger perspective. 


The beauty of studying the entirety of the function s(x), for X € 
ac, f € C?(G), x € G, which we shall explore in this paper, is that it 
could be viewed as a transformation in six (6) different ways; As 

(1.) z+ ky(A) = 8,¢(x), for any f € C?(G) 
and 
(2.) x +> ko(f) := sy, ¢(), for any  € ag, 


(from where the Plancherel formula for the space of functions x ++ ko(f) has 
recently been computed in [8.]) both of which are maps on G; or as 


(3.) f+ (A) := s(x), for any rz EG 
(which, at « =e, led Harish-Chandra to the consideration of f ++ (Hf)(A) : 
cf. [9.]) and 
(4.) f ++ lo(x) := 8),¢(x), for any » € ag, 
' both of which are maps on C?(G); or as 
(5.) AH mi(f) = 8),¢(“), for any r €G 


and 
(6.) A+ Mo(x) = 8y,¢(x), for any f € C?(G), 


both of which are maps on ag. Hence the function x ++ s),¢(x) may rightly be 
called an hyper-function on G whose major contribution to harmonic analysis 
would be to absorb other known functions of the subject and put their re- 
sults in proper perspectives, as we shall establish here for the Harish-Chandra 
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transform and Trombi-Varadarajan Theorem. 


In order to know the image of the spherical Fourier transform when re- 
stricted to C?(G//K) we need the following spaces that are central to the 
statement of the well-known result of Trombi and Varadarajan [9.]. Let C, 
be the closed convex hull of the (finite) set {sp : s € tw} in a*, ie., 


= {SoMa 20, yo =, nen} 
i=l i=1 
where we recall that, for every H € a, (sp)(H) = 4 Dyeqt my: X(87*H). 


Now for each € > 0 set §* = a* + icC,. Each F* is convex in ag and 


int(3) = (J ¥ 
O<e'<e 
((9.], Lemma (3.2.2)). Let us define Z(3°) = S(a*) and, for each € > 0, let 
2 (8°) be the space of all C-valued functions ® such that (i.) ® is defined and 
holomorphic on iné(¥*), and (it.) for each holomorphic differential operator 
D with polynomial coefficients we have supjp;ge) |D®| < 00. 


The space Z2(§*) is converted to a Fréchet algebra by equipping it with 
the topology generated by the collection, || - \|z(ge), Of seminorms given by 
| Dll zcgey = SUPsne(gey |DO|. It is known that D® above extends to a continu- 
ous function on all of §* ((9.], pp. 278 — 279). An appropriate subalgebra of 
2(§*) for our purpose is the closed subalgebra Z(§*) consisting of tv-invariant 
elements of Z(§*), « > 0. The following (known as the Trombs- Varadarajan 
Theorem) is the major result of [9.] : Let 0 < p < 2 and set € = (2/p) — 1. . 
Then the spherical Fourier transform f ++ f is a linear topological algebra 
isomorphism of C?(G//K) onto Z(§*). That is, the topological algebra 2(3*) 
is an isomorphic copy or a realization of C?(G//K). 


In order to find other isomorphic copies or realizations of C?(G//K) under 
the more inclusive general transformation map 
fOk(A) := 8) 7(2), for any x € G, 


we shall now introduce a more general algebra, Zq(¥*), of C—valued func- 
tions on int(¥*) x G which, when restricted to int(%*) x exp(Np), coincides 
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with 2(§*). The form of this new algebra is suggested by Theorem 3.5. Set 
ZG(8°) = S(a*) x G and let Z¢(¥), € > 0, be the collection of all C—valued 
functions V ((A,xz) +H V(A,z), V (A, 2) € int(F*) x G) such that 


(i.) U is holomorphic in the variable , analytic in z and spherical on G; 


(i2.) sUPinsge) [Di1V| < 00 and supg|VDz| < 00, for every holomorphic 
differential operator D, with polynomial coefficients and every left-invariant 
differential operator Dz on G and 


(iid.) the restrictien of V to int(¥*) x {e} (or to int(F*) x exp(No(A*)), 
for some zero neighbourhood No(A*) in g, as will later be seen in Theo- 
rem 3.5) is (a non-zero constant multiple of) the Harish-Chandra transform, 


(Hf)(A) = fF. 


It may be shown, in exact manner as for Z(¥*) above, that the space 
Za(8*) is converted to a Fréchet algebra by equipping it with the topology 
generated by the collection, || - || Zq(3¢), of seminorms given by 


Wllzewge) = sup |DVDzI. 
int($*)xG 


An appropriate subalgebra of Z¢(%*) for our purpose is the closed subalge- 
bra Z¢(§*) consisting of tv-invariant elements of Z¢(*), « > 0. By the time 
Theorem 3.5 is established it will be clear that 2;,.;(%*) ~ Z(§*), for every 


x in some zero neighbourhood No(A*) in g. In particular, 2;.}(%*) ~ Z(§*). 


3 Series Analysis of Spherical Convolutions 


Let f € C(G) and A € ag, we recall from [8.] the definition of spherical 
convolutions, s),f, on G corresponding to the pair (A, f) as 


sy p() = (f* yx), 2EG. 


We already know that s),¢(e) = (Hf)(A), where e is the identity element of 
G and A € ia*. This relation between a function on G at the identity element 
and another function on ia* suggests we study the full contribution of the 
Harish-Chandra transforms, (Hf)(A), of f to the properties of x ++ s),¢(z) 
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and to seek other functions on ia* which have not been known in the har- 
monic analysis of G, but still contribute to a deeper understanding of the 
structure of G. 


In order to explore the nature of this idea we consider opening up the 
spherical convolutions x ++ s),(x) via its Taylor’s series expansion. 


Lemma 3.1. Let No be a neighbourhood of origin in g andt be sufficiently 
small in R (sayO <t <1). Then 


oo 


ern 
8y,¢(zexptX) = > alt $y,f)(z); 


n=0 


where for every X € No we set [X"s),5|(x) = f-5y,f(Z exXPUX)),_o 
Proof. The proof follows from a direct application of Taylor’s series ex- 
pansion, [5.], p. 105. O 


At x =e and ¢ = 1 the formula in the Lemma becomes 


sy p(x X) = > Ray g](e) = saple) +) plan C6) 
n=0 n=1 


= (HAO) + LK sysH(O), X € No. 
n=1 


This observation leads quickly to the following result which gives the exact 
contribution of the Harish-Chandra transforms to the study of spherical con- 
volutions. 


Lemma 3.2. The Harish-Chandra transforms, \++ (Hf)(A), f € C(G), 
is the constant term in the (Taylor’s) series expansion of spherical convolu- 
tions, z+ 8,¢(x) around x =e, for everyX €a*. O 


It may be deduced, from the expansion leading to the proof Lemma 3.2, 
that the only time the remaining terms in s),s(exp X), after the (non-zero) 
constant term (Hf)(X), could vanish is when the differential operator X = 0. 
That is, when X = 0. It therefore follows that the well-known (Harish- 
Chandra) harmonic analysis on G ((1.], [2.], [9.] and [11.]) has always been 
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that of the consideration of the map X ++ s), ;(expX) at only X = 0, which 
is the origin of g or which corresponds to the identity point of exp(g). Hence, 
since the constant term, (H/f)(A), of 8), s(exp.X) corresponds indeed to the 
consideration of the constant term in the asymptotic expansion of (zonal) 
spherical functions, ~), it also follows that other terms in the expansion of 
py may be needed to completely understand f ++ s), ;(z). 


The expression for s), ¢(expX) therefore suggests that a full harmonic 
analysis of G may be attained from a close study of the remaining contribu- 
tions of the transform of f given as 


A SLR", se), 


for all X € No, ne NU {0}, « € G, f € C(G) and sufficiently small values 
of ¢, in the same manner that its constant term, 


A— (Hf)(A) 


had been considered. 


However before considering the transformational properties of spherical 
convolutions we note the following lemmas which lead to a more inclusive 
view of the Trombi-Varadarajan Theorem and prepares the ground for its 
generalization. 


Lemma 3.3. Let No be a neighbourhood of origin ing, A € ag and t be 
sufficiently small in R (sayO <t <1). Then 


8y,s(wexptX) = ey 


=o 


for every X EN, EG, f €C(G). 
Proof. We note here that 


= \(Diaal 8y,f(2), 


em 


(Xs, ,](@) = Fon (COP uX)oaa = SUS #Pa)(TOXDUX) oc 


= (f #29) (COD UX) ng = WLIO) (F # Ps)(CORP UX) ao 
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Hence 


(55 /](a) = o( 


du” 


Yio Fer) (#2) = 1G) Mc'sry()- 0 


The particular case of setting x = e and ¢t = 1 in Lemma 3.3 introduces 
the Harish-Chandra transforms, (Hf)(A), into the analysis of this series, 
proving the following. 


Lemma 3.4. Let No be a neighbourhood of origin in g, f € C(G) and 
A € a*. Then the spherical convolution function, © ++ s),f(x) is a non-zero 
constant multiple of the Harish-Chandra transforms, (Hf)(A), on exp(No). 
Proof. Set « = e and ¢t = 1 into Lemma 3.3 to have 


oo 


sag(emX) = 1D (Gin) Mel asl) = 1D “y(oo)(A)ucol HE), 


n=0 
with paak t(D) co = =1+ et a7 £-)(A) yoo] # B..0 


Let us denote the non-zero constant in Lemma 3.4 above by x. The fol- 
lowing theorem is a consequence of normalizing the spherical convolutions in 
Lemma 3.4. 


Theorem 3.5. (Trombi- Varadarajan Theorem for Spherical Con- 
volutions) Let 0 < p < 2, sete = (2/p) —1 and x € exp(No). Set 
fe(A) = 4s (x) for f € C?(G//K). Then the spherical convolution trans- 
forms f+ f, is a linear topological algebra isomorphism of C?(G //K) onto 
2(¥). O 


We recover the Trombi-Varadarajan Theorem for Harish-Chandra, trans- 
forms by setting « = e in Theorem 3.5. Indeed, Theorem 3.5 above says 
that every x € exp(Npo) (and not just z = e) gives a topological algebra 
isomorphism between C?(G//K) and 2(%*). However if x € G'\ exp(Np), for 
any neighborhood No of zero in g, Trombi-Varadarajan Theorem may not 
be appropriate and it may be necessary to seek a more general realization 
of C?(G//K) under the map f ++ i(A) := 5),f(x), for any c € G. Before 
considering another major result of this paper, giving the fine structure of 
spherical convolution functions, we state a result on the finiteness of a central 
integral usually used in the estimation of many other integrals of harmonic 
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analysis on semisimple Lie groups. 
To this end we define, for every x € G, the function z ++ d(x) as 
dle) = [ Sy-*2)(1 + o(y""s)) Pay 
We observe here that 
> d(e) = [ By" )(1 + o(y")) dy = iE 2’(y)(1 + o(y)) "dy, 


which is a constant whose proof of finiteness may be found in [11.], p. 281. 
This constant is crucial to all harmonic analysis of C(G) and, in particular, 
to the embedding of C(G) in L?(G). It is therefore important to understand 
the nature of d(x) for all x € G in order to employ it in a more inclusive har- 
monic analysis on G. We consider the nature of this integral in the following. 


Lemma 3.6. Let x € G. Then there exist r > 0 such that 


d(x) = cae + a(y~‘x))"dy < 00. 


Proof. We already know that =(y-!z) <1. Also 
1+o(y7'x) < (1+o(y))(1 + o(e)) = (1+ o(y))(1+ o(2)). 
It follows therefore that 


d(e) < i (14. o(yte)) dy < (1-+.0(a)) I (1+ o(y))dy. 


The last integral in the above inequality is finite if we embark on its compu- 
tation via the polar decomposition, G = K -cl(At)- K, of G. 0 


Theorem 3.7. Let No be a neighbourhood of origin in g where f ts a 
measurable function on G which satisfies the general strong inequality. The 
integral defining the spherical convolution function, x ++ sy,7(x), is absolutely 
and uniformly convergent for all x € exp(No),  € ia*. Moreover the trans- 
forms A+ s),¢(z) of f, with x € exp(Npo), ts a continuous function on ia*. 
Ifr > 0 is such that d(x) = [4 =(y~'z)(1+a0(y'z)) "dy < 00, z € G, then 


| sy,¢(2) |S d(x) + prap(f), 2 €G, A € ia". 


a6 de 


Proof. We recall that | y)(x) |< go(x) = E(x), x € G, A € ia*. Hence 


| (frpr)(2) |S a | f(y)ea(y*s) | dy S Maar (f) [ =*(y'a)(1+0(y*2)) "dy 


= d(x) - f1,1,r(f). Continuity follows from the use of the Lebesgue’s domi- 
nated convergence theorem. 


The following well-known result on the foundational properties of the 
Harish-Chandra transforms, \ + (Hf)(A), A € 2a*, now follows from the 
general outlook given by Theorem 3.7. 


Corollary 3.8. ({9.]) Let f be a measurable function on G which satisfies 
the strong inequality. The integral defining the Harish-Chandra transforms, 


(Hf)() = a f(a)or(w)de, 


is absolutely and uniformly convergent for all A € 1a* and is continuous on 
ia*. Ifr > 0 is such that d = [,=*(y)(1 + o(y)) "dy < 00, then 


(HF)(A) |< dere (f), A € ia”. 


Proof. Set X = 0in Theorem 3.7 to have the first results. The inequality 
follows if we set « = e and observe that d(e) = [,2?(y~')(1+o(y1)) "dy = 
d. 0 


We now consider the image of C?(G//K) under the full spherical convo- 
lution map, f +> 1;(A) := s),;(x), for any x € G, In order to discuss this we 
have two options. One of the options is to introduce wave-packet that will 
still have its domain as Z(%*) while using an appropriate Plancherel measure 
on §*. This option has been explored in [8.], p. 34, where the L? Plancherel 
measure, d¢z,, on §* for the spherical convolution function (when viewed as 
a function on G) was defined to absorb the group variable, z. The results 
therein suggest that the image of C?(G//K) under the full spherical convo- 
lution map is indeed possible. 
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The second option is to retain the spherical Bochner measure, dX, on (a 
subset of) © and define the wave-packet as a map on the Fréchet algebra 
Z¢(§*). This will reflect the nature of the full spherical convolution map as 
a transform of members of C?(G//K) whose arguments are (generally) taken 
from int($*) x G (and not just from int(§*) as in the first option). This is 
the option we shall explore in the present paper. 


To this end recall the Fréchet algebra Z¢(3*), Ve > 0, let U € Zc($°*) 
and set 


No(At) = Non At, 


where No is a zero neighbourhood in g. It is clear that No(A*) is also a zero 
neighbourhood in g and that U = W(A, x), for all (A, x) € int(¥*) x G. It fol- 
lows, from Theorem 3.5, that 2;.}(%*) ~ Z(¥), for every x € exp(No(At)). 
We then have the following. 


Lemma 3.9. For every x € exp(No(A*)) and U € ZG(8*), we have that 
UA, z) = (A), for some ® € Z(F*). 


We now employ these remarks to define a map from Z¢(§*) to C?(G//K) 
as follows. Let a € Zg(§*) and » + c(A) be the Harish-Chandra c—function 


defined on $; := ia*. We associate to every a € Z@(¥*) the function y, on G 
defined as 


pa(e) =| tw |n# [ a(—h2)eilaed=) “e010. @ EG. 


It should be noted here that 


Qa(z) =| to | [ a(—A, x)p_y(x)e(—A)~*e(A) “dA 


=| 10 | [ a(A, ©) py(a)e(d)“te(—A)-*d(—2) 


=| to |~! i a(A, z)py(x)e(A)~1e(—A)~"dd, 
$1 
which is due to the invarianve of dA, and that 


Palki zk) = va(z), 
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Va€G, ky, ko € K, being a property inherited from a and 9). 


The (extra) requirement of being spherical on G placed on members of 
Ze(8°) may at first be seen as a restriction, when compared to the require- 
ments on members of 2(¥*). It however turns out that this extra require- 
ment is what is needed to assure us of the generalization of the classical 
wave-packets (of Trombi-Varadarajan) on G to all of t +> (x). This is 
established as follows. 


Lemma 3.10. Leta € Z¢(%*) and No(A*) be as defined above. Then, for 
every z € exp(No(A*)), the map x ++ Ya(z) ts the classical wave-packet of G. 


Proof. We observe that, with exptH € exp(No(A*)), 
a(A, ©) = a(A, ky exptHke) = a(A, exptH) = ®()), 


for some © € 2(§*). Here we have employed the spherical property of a on 
G in the second equality and Lemma 3.9 in the third equality. 0 


The above Lemma shows that the definition and properties of the map 
LH a(x), x € G, is consistent with the relationship (in Lemma 3.4) existing 
between spherical convolutions, s(x) and the Harish-Chandra transfroms, 
(Hf)(A). Hence in order to extend Trombi-Varadarajan Theorem (which 
gives the image of the algebra C?(G//K) under f + (Hf)(A)) to alla eG 
(under the spherical convolution tranform), it will be necessary to show that 
Z +> %a(x) is the wave-packet of f ++ s(x) for all « € G. According to 
Lemma 3.10, this needs only be done for those « = ky exptHk, in G with 
exptH ¢ exp(No(A*)), for any neighbourhood, No, of zero in g. We however 
give a self-contained discussion of these results, the first of which is given 
below. 


Theorem 3.11. y, € C?(G//K) for every a € Ze(8*). 


In order to establish this Theorem we prove some lemmas which give 
appropriate background for it. Indeed we derive an appropriate bound for 
| ~a(h;u) |, where u € L(g) and A is well-chosen, and the appropriate col- 
lection of seminorms are also in place. 
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4 Algebras of Spherical Convolutions 


We now consider the various algebras of spherical convolutions that have 
emanated in the course of this research and their relationship with the Harish- 
Chandra Schwartz algebra, C(G), on G as well as its distinguished commu- 
tative subalgebra, C(G//K), of (elementary) spherical functions. 


Define C,(G) = {sy : f € C(G)} and set Cyo(G) = {s),,,}, for all 
A € ag. It is clear that U,. af C,(G) is contained in C(G). We may therefore 
topologize near C,(G) by giving it the relative topology from the topology 
defined on C(G) by the seminorms, [io,b,r- 


Lemma 4.1. The inclusions 


[U Cxo(@)]¢ e(@//K) c [UJ e(@)] < e(@) 


Eng, Eng, 


are all proper. 0 


Theorem 4.2. Uc: C,(G) is a closed subalgebra of C(G). 
Proof. We recall that [aoe(f * Yr) S Cli,byrtro(f) * Me,tyr(Yx), where 
c:= [,=*(x)(1+0(x))-da < oo for some ro = 0. However 


Ma,tir (Pr) = supllea(1; a; a)| -E(x)~*(1 + o(2))"] 


=| 1(@)(A) | -supllen(@)| -2(@)1 + o(@))] 
<M | 4(a)(2) | -sup[E(e)*(1 + o(@))"] < 00 


(since ~) is bounded for all \ € '). 
Hence Mapr(f *x) < 00, VAE FF. O 
It may be recalled that members of C(G) are exactly those functions on 


G whose left and right derivatives satisfy the strong inequality. In the light 
of this observation we define C(G) as exactly those functions on G whose 
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left and right derivatives satisfy the general strong inequality, for each « € G. 
Explicitly we set C™(G) as 


COG) ={f: GHC: supll F(a; y;b)| -E(y 2) "(1 + o(y*x))"] < co}, 
x € G. A collection of seminorms on each of C)(G) may be given by 
uo) (f) = supllf (a; y;b)| Bly *a) "(1 + o(y*2))"]. 
It is however clear that C)(G) = C(G), so that C(G) C Useg C/G). 


Theorem 4.3. The natural inclusion UsegC™(G) C L?(G) has a dense 
tmage. 

Proof. It is known that the natural inclusion of C(G) in £7(G) has a 
dense image, [{1.]. The result therefore follows if we recall that, as sets of 
functions, 

c(a@)c LJc@@) < P@), 
«eG 
where the second inclusion holds from the fact that d(x) < oo, cé€G. O 
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The relativistic quantum equation is proposed for the complex 
wave function, which has the meaning of a probability amplitude. 
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1 Introduction 


Fundamental issues related to the completeness of the description of reality 
in quantum mechanics were discussed long ago in well-known works [1,2] 
and then in many others. These discussions continue in present time. The 
Schrédinger equation for the complex wave function ~, which, according 
to Born’s interpretation, has the meaning of a probability amplitude, is 
nonrelativistic. This means that the theory allows propagation of signals at 
arbitrarily high speeds. Perhaps, some of the discussed difficulties and para- 
doxes of quantum mechanics are related to this circumstance. It should also 
be noted that quantum theory is used in the description and interpretation 
of experiments with photons, which are relativistic objects. So far there 
is no generally accepted relativistic equation for the field, which could be 
interpreted as a probability amplitude and would allow a Born probabilistic 
interpretation. The probabilistic interpretation of the complex field within 
the framework of Dirac theory was discussed in works [3,4]. 

In this paper the relativistic equation is proposed for the complex scalar 
field, allowing its physical interpretation as a probability amplitude. The 
theory is formulated within the framework of the Lagrangian formalism. 
The spreading of a wave packet in an unlimited space is considered. The 
relativistic correction has been found in the theory of a harmonic oscilla- 
tor, leading to a violation of the equidistance of levels. The issue of the 
completeness of the quantum description is not addressed in the paper. 


2 Relativistic equation for the complex scalar 
field 


The recipe for the transition from the classical to the quantum mechanical 
equation consists, as is known [5], in replacing momentum and energy in 
the classical Hamiltonian H = p?/2m with the operators p + —ifAV and 
H — ih0/0dt, which act on the complex function 7. In the relativistic 
case, the relationship between energy and momentum is given by formula 
H = V/m?c! + cp?, where m is the particle mass, c is the speed of light. 
Using the indicated substitution, we arrive at the relativistic equation 


nee =ovmc —HeA 7, (1) 


='63%s 


where A is the Laplace operator, o = +1. This equation is inconvenient 
for a number of reasons. Firstly, the time and space coordinates enter here 
unequally, so that it does not have an explicitly covariant form. Secondly, 
because in Eq. (1) the Laplace operator appears under the sign of the root, 
and, therefore, this equation is nonlocal. In order to get rid of the noted 
formal shortcomings, the transition to the quantum mechanical equation is 
usually carried out in the expression for the square of energy H* = m?c4 + 
c’p?, which leads to the Klein-Gordon-Fock equation 


160%" 43 
(a-355-2 ‘\y=0, (2) 


where A = hi me is the Compton length. In the following we will also use 4- 
dimensional notation, in which an arbitrary 4-vector is B = B, = (B, By = 
iBo), in particular ¢ = a, = (x, L4 = 129 = ict). The scalar product 
of two vectors is written as AB = A,B, = AB+ A,B, = AB — ApBo. 
A summation from 1 to 4 is implied over repeating Greek indices, and a 
summation from 1 to 3 — over repeating Latin indices. In this notation 
equation (2) can be written in the explicitly covariant form 


& _ x) Ho = 0. (3) 


From this equation there follows the continuity equation 


Op os 
where 
= A Oy On* ° : R * _ * 
pmig(y ca +), ie ee) 


Here the quantity p is not positively definite, and therefore cannot be inter- 
preted as a probability density. In addition, the transition to a higher order 
equation leads to the emergence of new solutions. In connection with this, 
the problem arises of constructing the quantum relativistic equation for the 
complex field, which would have the meaning of a probability amplitude. 
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3 Relativistic quantum equation for a prob- 
ability amplitude 


In spite of the shortcomings noted above, let us consider equation (1) as the 
equation for the probability amplitude ~(x). The advantage of this equa- 
tion, which to a great extent compensates for the difficulties noted above, 
consists in the possibility of interpreting the complex field as a probability 
amplitude. In this equation the sign in front of the root can be chosen 
arbitrarily, since this will not affect the physical results, but only affect 
the form of the time dependence of the wave function for stationary states 
wy ~ exp(— io Et/h). In order for equation (1) to have the generally ac- 
cepted form of the Schrédinger equation in the nonrelativistic limit, we will 
assume in (1) o = +1. The equation 


ine = V/mict — RA (6) 
does not have an explicitly covariant form, but leads to the continuity equa- 
tion for the probability density |9|?: 


Alp? — 7 
ava +V-j=0, (7) 


where 5 
V jai (wi HA pov -¥7A vy, (8) 


and j has the meaning of the probability flux density. The continuity equa- 
tion (7) for the probability density can be written in 4-dimensional form 


Oju(@) _ 
Ox, LS 0, (9) 
where j, = (j, tcl?) is a 4-vector of flux density, so that it has the same 
form in any inertial reference system. Thus, the equation for the probability 
amplitude (6), written in an arbitrary reference system, leads to the Lorentz 
covariant form of the probability conservation law (9). Since under the 
Lorentz transformation x, > i, = Aw%y, Where QyyQyy" = Oyy', the 4-vector 
is transformed according to the law j/,(z') = aujv(x), then the probability 
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density and the probability flux density in different systems are connected 
by the relations 


icly'(2')|? = iclep(x)|? + Gamjm(Z); 
Hula!) = ana telb(2)|? + dimijm(2)- 


(10) 


Let us represent the complex field as the sum of the real and imaginary 
parts 


1 / oon 
w(x) = wa (x) + iy"(z)]. (11) 


Then equation (6) turns out to be equivalent to a system of two equations 
for the real fields 7)’ and w”: 


he! = JV mit — Re by", nb = Vm A-ROD (12) 


Here and henceforth we also use the notation o = 0w/dt. One of the 
functions, for example w”, can be eliminated, and then we arrive at the 
explicitly covariant Klein-Gordon-Fock equation for the real part 


Wy — 7? Ad! +2 y! = 0. (13) 


Along with this, the evolution of the imaginary part is determined by the 
second equation (12). Similarly, the Klein-Gordon-Fock equation can be 
obtained for the imaginary part of the function (11). 

Let us consider the meaning of the expression 1 — 47A 7, containing 


the square root of the Laplace operator. For 2? < 1 the following expansions 
are valid [6] 


Vi£a=1-)0(F)"an2”, 
n=l (14) 
= 20 = 3)! Be gai oe 
ee ah ek gee gh ge 
(Lata)? =14+ o(F)"bn2”, 
n=l (15) 
» = (Ma)! ee ee ee: 
weer (Ombre vee gee eg os age 
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where (2n)!] = 2-4-6...2n = 2"n!, (2n+ 1)! =1-3-5...(2n+1) = 
Qn 


Te (n+3 / 2), With allowance for these representations, by the expression 
with the square root from the operator we will understand the following 


expansion 
V1—-VA p= (1 ->> anda") (16) 
n=1 


This expansion is valid under the condition |A7A?| < 1, which is assumed 
to be satisfied in the following. This means that one considers the states in 
which the wave function changes over distances greater than the Compton 
length. The opposite ultrarelativistic case can be considered in a similar 
way. It should be noted, however, that it is of rather methodological interest, 
since in this limit the possibility of transformation of particles should be 
taken into account, which can only be done by passing to the representation 
of secondary quantization. 

Taking (16) into account, we obtain a formula for the probability flux 
density. The divergence (8) can be written in the form 


s me = 2n (ph* AT Ny ly* 
Vi j= ime Dann (Wt and — pA"y'). (17) 


n=1 


From here it follows that the probability flux density has the form 


[oe) 
j= 05, 
n=1 


J) = ica Xt Y > [(Aemtys) (VA) — (ATH) (VAr*y") | 


e=1 
(18) 
Everywhere it is assumed that the Laplacian to the zeroth power is equal 
to unity A° = 1. 
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4 Lagrangian formalism 


Equation (6) can be obtained within the Lagrange formalism, if we choose 
the Lagrangian in the form 


L= ©" (itd - meV —¥A ) - © (init +méV1—A yp"). 
(19) 
Lagrangian (19) is a function of variables », y* , , 4" and derivatives with 
respect to spatial coordinates Ay, A?y,... A?w,... Ay*, A?y*,... A?y*,. 
In this case, the Euler-Lagrange ecuntions Wave the form 


OL OOL OL gob 7 
Bp as aa 70s ane ee (20) 
and a similar equation with the replacement 7) > 7*. A substitution of the 
Lagrangian (19) into (20) leads to Eq. (6). 
In the absence of time-dependent external fields, from the requirement 
of invariance of the Lagrangian under the time shift t > t+ t and w(t) > 
w'(t +t), there follows the continuity equation for the energy density 


OH 
= 21 
OE — +Vsjz =), (21) 


where the energy density 


eOLE NE Te J a rere 
=y— sith 5 ah La ly v1 MA b+ ov XA y"], 
(22) 
The energy flux density is given by the formulas 
n=1 
IP = 5° [(arey)(Wartz,) — (Wary) (At) + 3) 


a=l 


+a") (VASILE) — (VAT) (Ae Zs)]. 
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Here we used the notation 


OL y* 6L w 
L,= =meaAr"—, DL =——— = me'a,r""— 
An mea 5 Ey Amp MC" nA 5 be 
AM he = metagrm ae, ASL. = meta,xr =. 


Taking (24) into account, we have the final expression for the energy flux 
density 


j2 = nc anX?” S [(a"-*4) (VAc-1y*) _ (VA"-%p) (A*1y*) ae 


a=l1 


+(A-ey)") (VAc-1p) ae (VA"-=7h*) (ax-*y)]. 


(25) 

The condition of invariance of the Lagrangian with respect to spatial 

translations x + x’ = x + Xo leads to the continuity equation for the 
momentum density: 

T;, + Ve Oink = 0, (26) 


where the momentum density 2; and the flux density of ¢- component of 
momentum oj, are given by the formulas 


OL OL th 
~=—V *—— — V; —— Se *V; = V; *), 
T air vay 5 (b'Vid — vVi¥") (27) 
Ck = ce 
n=1 


of”) Si 3 [(a”*vi0) (V,AeL,) = (ViA"-*V 2h) (At) + (28) 


a=1 
+(A™ Vit) (Ve? ILs) — (Vad oVi") (ATLA). 
With allowance for (24), the last formula takes the form 


o® = —sncan x x 


x S [(a-°V.0) (Vide tp") — (ViA™ Vip) (At Ty") + (29) 


a=1 


+(A™° Vit") (Vado ty) — (Ved *Vid") (At p)], 
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Lagrangian (19) is also invariant under the phase transformation 


V(x) + $(e) = v(a)e™, (30) 
where the parameter a does not depend on coordinates and time. A conse- 
quence of this phase symmetry is the probability conservation law 


Ob _ yy Ob. 
(¥55 Fi ae) 
+Vi DD [(Vid ey) ATL, — (AN) VAT LA (31) 


n=1 a=1 
= (VA eap*) Ae IL * ie (Arey*)ViAe zs] acify. 


This equation, with taking into account formulas (24), coincides with the 
continuity equation (7) for the probability density which is obtained directly 
from equation (6). 


5 Interaction with the electromagnetic field 


The interaction with the electromagnetic field is enabled using the well- 
known derivative substitution 
6) 6) .e 
dn, card Oa, + te 
where e = +le|, A, = (A, Ag = i®) is a 4-vector potential of the elec- 
tromagnetic field. In three-dimensional notation, the substitution (32) is 
equivalent to the following substitutions: 


An; (32) 


e 6) 6) 

ee pe ck ge rae 
VaIV+e A, is, + the +e ; 
e 2e e\? 9 (33) 


As a result, we obtain the quantum relativistic equation for the probability 
amplitude in the electromagnetic field: 


ince —U(2)p = Vmict — Re? Ay ¥, (34) 


where U(x) = —e®(z) is the potential energy in an external scalar field. 
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6 Schrodinger equation with the relativistic 
correction 


Let us write the Schrodinger equation with the scalar potential, taking into 
account the main relativistic correction. In this approximation 


1 
V1~VAys (1 ~ 5X A - 5a") yp. (35) 


Then equation (34) will take the form 


Oe af, 132, _ 1yan2 
ine = Ua) +m 5A 304 w. (36) 


In order to proceed to the nonrelativistic limit, it is convenient to pass 
to the wave function x(x), which differs from the original function by a 
time-dependent phase factor 


mo* 
U(a) = x(2)er*. (37) 
Then equation (36) will take the form 


Ox i? iit 2 
ii = —-— —- — A’. 38 
Ot 2m Ax + U(x) x Bm * ee) 
Here the last term is the relativistic correction. The probability flux density 


(18) j = j@+j® is the sum of the nonrelativistic part j and the relativistic 
correction j): 


ih : 
jf = -—— (x*Vx - xVx"), 
je ant (39) 
j? = — a (x*VAx — xVAx* + Ax* Vx — AxVx*). 


The energy density (22) and the energy flux density jz = j2 + ji? (25) 
in the considered approximation are given by the formulas 


ee x ‘ Z 
H = me|x|? — ran Ax + xAx*) - a: A’y + xA*x*), (40) 
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is) = 


WP ok ‘ * *kO. a h 2 * * 
re 6 Vx +XVX" = XVX — XVX" + tg me (x*Vx — xVx"*), 


iP = Fm [(vx7) (av) + (Ax) (Wx) + (VAY) X+xX'(VAx)- 
—x* (VAX) - (VAx')x — (Ax) (Vx) - (Vx) (Ax)]- 


2 
~iF me?x?|(Vx") (Ax) — (Ax’) (Vx) + (VAxt)x— 22 (VAx)]. 
(41) 
As we see, in the continuity equation for energy, even in the nonrelativistic 
limit, it is necessary to take into account the rest energy of the particle. 
The momentum density (27) is expressed through the probability flux 
density in the nonrelativistic approximation by the relation 7 = mj. The 
flux density of 7-component of momentum o;, is given by the formulas 


1 2 
Oi = of a cae 


hi? 
Pil = hee, Fa (ViVi x) +x(ViVex*) — Vix* Vax — VixVex'], 


x Ix" (ViVe Ax) +x(ViViAx*) + Ax* (ViVi Ax) + Ax(ViV Ax’) — 


—(ViAx") Vix — (ViAx*) Vix — Vix" (Vidx) — Vex" (V:Ax)| 
(42) 
In equation (38) the last term should be considered as a small perturbation. 


7 Evolution of the wave function in an un- 
limited space 

Let us consider the evolution of the wave function in an unlimited medium. 

Let us move on to the new wave function y(z) (37), such that (2) = 


x(x) exp(—it/tc), where t¢ = X/c = h/mc*. Note that for an electron 
To & 1.3-107?4s. In this case equation (6) will take the form 


itcx = (vi =VAq 1)x. (43) 


Ber pie 


The equations for the real and imaginary parts of the wave function y = 
(x! + ix”) / V2 can be written in the form 


Tox! = (v1 — PA - 1)", toxl = -(v1 =A = 1)y, (44) 


There holds the normalization condition 


= [a x(x, t) + x" (x, t)| =i, (45) 


Let us consider the evolution of the wave function in an unlimited space, 
provided that at the initial moment ¢ = 0 there are given the functions 
x'(x,0) and y”(x,0), normalized by the condition (45) . From equations 
(44), it follows a second order equation in time for the function y’: 


TAX! + (v1 —VvA— 1) x =i (46) 


We expand the required functions into the Fourier integral 


1 ikx 1 ikx 
x'(x, t) = (anye [aoe dk, x"(x,t) = aaa | uW« dk. 
(47) 


Substituting these expansions into (46) and the second equation (44), we 
obtain 


o 


x +wExe=0, Xx" + wExk = 0; (48) 
where 
we = (V+ 0H? - 1) Ee (49) 
To 
In the nonrelativistic limit Ak < 1 
Ak? Tike 
wa SE Me (50) 


The group velocity, which determines the speed of propagation of a wave- 
particle 


V ae ee eee 51 
oO dk Jy RR” ey 


cannot exceed the speed of light. A calculation of the group velocity within 
the framework of ordinary nonrelativistic quantum mechanics by formula 
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(50) gives v, = Ak/m, so that with an increase of k or a decrease of the de 
Broglie wavelength the speed of signal propagation can be arbitrarily large. 
Solutions of equations (48) have the form 


Xi,(t) = Ape + Bye? xu(t) = —i Ape" + iBpe™s?, (52) 


From the conditions for reality of the functions x’, x” it follows vy, = x’%*, 


and x4 = x/%, so that B, = A*,. The normalization condition (45) gives 


[ae (uP +A) =, (53) 
The coefficients in (52) can be expressed through the initial conditions 


Ae = $(x4(0) +140), AL. = SOKO —ixk)), (54) 


so that ; ; 
Xn(t) = X%(0) cos wpt + X4,(0) sin wyt, (55) 
55 
Xe(t) = —Xj,(0) sin wet + x4 (0) cos wet. 


Taking into account the relations 


1 —ikx 1 -i 
X4(0) = aan | ae *y'(x, 0), x%(0) = aan | dee ""(x, 0), 
(56) 


we find solutions which determine the real and imaginary parts of the wave 
function at an arbitrary moment of time from their values at the initial 
moment of time 


x’ (x,t) = [x [Jo(x — x’, t)x’(x',0) + Js(x — x’, t)x"(x’, 0)], 


(57) 
x" (x, t) = [et — Is(x — x’, t)x!(x',0) + Jo(x — x’, t)x"(*', 0], 
where ; 
Jo(x — x,t) = an ju cos wt eRO*) 
(58) 


Jg(x —x’,t) = on [us sin wyt eb*O-*) | 
1 


2942 


Using the expansion of a plane wave in terms of spherical functions and 
integrating over angles, we obtain 


1 = 
Jo(x — x’, t) = Jo(|x — x'|,t) = =a | dk k? cos(wyt) Jo (K|x — x'l), 
0 


1 ri ; ‘ 
Js(x— x,t) = Is(be- x11) = 55 i dle k? sin wt) Jo (lx — x1), 
(59) 


where jo(x) = sinz/z is the spherical Bessel function. 
Let us express the sought functions at the initial moment ¢ = 0 in terms 
of the modulus and phase 


x (x, 0) = V2 |4h(x, 0)| cos (x, 0), x(x, 0) = V2 |w(x, 0)| sin O(x, Meg 


As the initial condition we choose the Gaussian distribution for the proba- 
bility density : 

[(7, 0)? = aa (61) 
where r = |x|, and the parameter co determines the width of the wave packet, 
so that as o decreases the distribution (61), satisfying the normalization 
condition 47 f5° |y(r,0)|?r?dr = 1, approaches the delta function. For the 
initial phase in (60) we choose 6(x, 0) = 0, so that for the real and imaginary 
parts we have 


a er 1a"), x" (r, 0) = 0. (62) 


x'(r, 0) = V2 |(7,0)| = 


In the calculation it is helpful to use the well-known integral [6] 


oO v—1/2 2 
= 242 v+3/2 - = vis b _ 2» 
| ert ¥ jv—2(bt) dt = oan e€ 407, (63) 
As a result, we obtain a solution in the form 


[e,0) 2 2 
x'(r,t) = ae =o A i dk k? cos(wyt) Jo(kr) e7 "=, 
0 


x"(r,t) = ee 2 / dk: k sin (wt) jo(kr) cE. 
ss 0 


(64) 


vi 


a: 75:= 


Figure 1: The spatial distribution of the probability density |p|? = 
n93\|? for o/X = 1 at the moments of time tr = t/te : 1-0; 2- 
0.5; 9- 1.3; 4- 2.0. 


The spatial distribution of the probability density |p|? = 7¥%o3|y|? = 
n¥3|y|? at some moments of time is shown in Figure 1. As we see, a 
particle localized at the initial moment of time at the origin of coordinate 
system is gradually with equal probability spreading throughout space. 


8 Relativistic correction in the oscillator 
theory 


The Schrédinger equation for a one-dimensional along the z-axis quantum 
oscillator with account of the main relativistic correction (38) has the form 


OX i? Oy mus? it omy 


— = -— 5, 65 
OE 2m Ox? 2 ~~ 832 dat ee) 

In the stationary case, when (2, t) = y(x) exp(—iEFt/h), we have 
_ Way mura? ht dty (66) 


“Om de® "2° B82 det 
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It is convenient to introduce a dimensionless coordinate = (/mw/ h)a, 
and in result the equation takes the form 


ep = —p") + Py — yp), (67) 
where differentiation occurs with respect to £, ¢ = 2E/hw, and 


fw 

=-——,. 68 
Y= ime (68) 
In order to estimate this dimensionless parameter, we take iw ~ 1K & 

10~!% erg, m + 10-78 g. As a result, we get a very small value y ~ 107", 
The correction to n-th energy level is determined by the diagonal ele- 
ment of the perturbation operator taken over the unperturbed wave func- 

tions [5] 


” 
eA) = —(n| Fi |n). (69) 


The wave functions of an oscillator in the coordinate representation are 
expressed through Hermite polynomials [5] 


mu \ "4 a 
(sin) = (SP) ne Fa. (70) 


As a result, we find 
=~ 75 (n+1)(n+2). (71) 


Thus, the energy of the oscillator level with account of the main relativistic 
correction takes the form 


By = ta] (1~ $7) + 51-39) — Fam? ; (72) 


As we see, the relativistic effect leads to a decrease of the ground and excited 
levels, as well as to a violation of the equidistance of the spectrum. The 
effect of violation of the equidistance property, leading to a decrease in the 
distance between higher levels, can be used for the experimental detection 
of relativistic effects in quantum mechanics. 
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9 Conclusion 


The paper. proposes the relativistic generalization of the Schrodinger equa- 
tion for the complex function, which can be interpreted as a probability 
amplitude. Unlike the nonrelativistic equation, this equation allows prop- 
agation of signals only at a speed not exceeding the speed of light. The 
theory is also formulated within the Lagrangian formalism. The equations 
of conservation of the probability, energy and momentum of the complex 
field are obtained. The problem of evolution of a wave packet in an un- 
limited space is solved. The relativistic corrections to the energy levels of 
the harmonic oscillator are found. It is shown that relativistic effects lead 
to a shift and a violation of the level equidistance, which can be used for 
experimental detection of the relativistic effect in the probabilistic quantum 
theory. 
The author thanks A.A. Soroka for help in preparing the article. 
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The hypergeometric functions occur naturally in variety of problems 
in analysis, statistics, operational research, theoretical and mathematical 
physics and engineering sciences. In the physical problems where these 
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and integral representations of them and discuss several simply-structured 
single term approximation formulae: for the Hubbell type radiation inte- 
grals, generalized hypergeometric functions and Appell’s functions. We 
propose also a new approach to estimate the pFy-functions by relating 
them to three simple elementary functions. The notions generalized (mul- 
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1. Introduction 


The generalized hypergeometric functions are defined by the series 


({3],[17]) 


Q1,42,...,@ 
Ce ee ae ey eae 
. k (1) 
=a COE , where (a), =I'(a+k)/T(a) . 
k=0 (br) (B2)e-n(ba)e RE 


It is well known that the pFj-series converges for all finite z ifp<q, 
for |x| < 1 if p=q+1 and diverges for allxe 40 ifp>q+1. 


According to [16], we consider three separate classes of pFy-functions: 


e p<q. Every ofj-series is a Bessel function. Special cases include 
the trigonometric and hyperbolic sine and cosine functions. The gener- 
alized sine and cosine functions of order g and the so-called hyper- Bessel 
functions J; (2) .¥g(£) ([1]) are also functions of this class, beeing: oFy-1. 


We call all such functions generalized hypergeometric functions (g.h.f-s) of 
Bessel type. 


e p=q. The simplest examples are: ofp — the exponential func- 
tion and ,F, — the confluent hypergeometric series. Their special cases 
include the error functions, incomplete gamma functions, exponential and 
logarithmic integrals etc. Cases in which the series terminates include the 
Hermite, Laguerre and Bessel polynomials. All these special functions will 
be called generalized hypergeometric functions (g.h.f-s) of confluent type. 


e p=qt+l1. The ,Fo-function is the binomial series and 2F, is the 
well known Gauss hypergeometric function. Its special cases include loga- 
rithms, inverse trigonometric and hyperbolic functions, Legendre, Cheby- 
shev, Gegenbauer and Jacobi polynomials and functions, incomplete beta 
functions, complete elliptic integrals of first and second kind, etc. Consider 
them as a class of the Gauss type g.h.f-s. 


Exton [4] discussed a number of problems involving finite and infinite 
statistical distributions as well as a number of engineering problems which 
gave rise to integrals associated with hypergeometric functions of one and 
two variables. The enormous success of the theory of hypergeometric series 
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of one variable stimulated the development of a corresponding theory of 
two and more variables. Thus, we have the Appell’s series which are 
defined in the following form [3]: 


[o,e) / 
b b m n 
Cet eee (2)m+n(>)m(2 )n 2 a , max{|z|, ly|} <1: 


eral (c)min m! 
. (2) 
(2)m+n())m(b )n 2™ 


ee n 
/ / y 
F(a, bb 3c,¢;2,y) = Ss; " a 5 |z|+|y| <1; (3) 


m,n=0 (C)m(C)n m! 
ae s+ (a)m(@ )n()m(B Jn 2 y 
eee RE : 
F3(a,a ,b,b;c;2,y) = > is al »max{|z|, |y|} < 1; 
m,n=0 


(4) 
[,@) 
b gm y" 
Fi(aybjecs2,y) = Y) Amsnllmen PV Hale Vivi <1, (6) 


m,n=0 (c)m(e’)n m! 


; ! : 
where as usual, the denominator parameters c and c are neither zeros nor 
negative integers. 


For the Gauss hypergeometric function the so-called Euler integral rep- 
resentation is well known ([3],{17]): 


2Fi (a, 837; 0) = rariecs [e — uy? PtP (1 — eu)-%du . (6) 


0 


Appell gave some double integral representations of Euler type for the 
series F\, Fy and F3, for example ([3)): 


F,( Poe CR aoa ge 1Y) T(B)I(B P(y — BYE’ 4B) 


ica 
/ U / 
x [ [ew “11 — u)?~P-101 — v)1 —F 11 — ue — vy) 7% du dv , 
0 0 


(7) 
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provided Re(y) > Re(f) > 0, Re(7 ) > Re(/’) > 0. 
Another representation of the same type is ([3]): 
ry’) | fm 
POG Aran —a) J J 


x (1 — 0179-11 — va) ~% w= w)t 811 = ve - wy)? dv dw 
(8) 


t / 
F(a, 8,8 57,752,y) = 


- with Re(y) > Re(S) > 0, Re(7’) > Re(a) > 0. 

Several other integral representations have been found for the Appell’s 
functions by various authors, all of them like (6),(7),(8) involving fractional 
type integrals (single or double, see Section 5 and [16],{19]). Picard pointed 
out that the series F, can be represented also by a single Eulerian type 
integral: 


(7) 


F(a, 8,872.9) = aE ay 


1 
/ 
x pera Say = ux) F(1 —uy)~? du , Re(y) > Re(a) > 0, 

0 

(9) 
however there appears to be no single integral representation of this type 
for the Appell series Fy. Burchnall and Chaundy gave a double Eulerian 
integral for this function, similar to (7). 

In all these representations it is assumed that |z| and |y| are small 
enough to make both series and integrals convergent. 

There are many approximation formulas of different type, for the spe- 
cial functions of mathematical physics, and especially for generalized hy- 
pergeometric functions and integrals involving them. Polynomial, rational 
and Pade approximations for (1) are well known and often used (see e.g. 
Luke [17]). However, we would like to survey some single term approxima- 
tions for hypergeometric functions and integrals, obtained recently by the 
authors and allowing fast and simple advanced evaluation of their values 
for practical use. 

In Section 5 we propose also another approach to estimate the gener- 
alized hypergeometric functions by relating them to three simple elemen- 
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tary functions. The notions generalized (multiple) fractional integrals and 
derivatives ({16]) are used. 


2. Single term approximations of radiation integrals 


e Hubbell rectangular source integral. Hubbell et al. [9] have 
obtained a series expansion for the calculation of the radiation field gener- 
ated by a plane isotropic rectangular source (plaque), in which the leading 
term is the integral 


b 
a a dz 
h(a, 6) = an [acts =) VEEL a> 0, b > 0. (10) 
0 


The above integral, called Hubbell zntegral, has found important appli- 
cations in metrology situations where mechanical contact is not allowed. 
Typical examples can be found in a variety of medical, agricultural and 
industrial applications which make use of radiometric gauging and pro- 
cess control. In fact, these and other metrology devices are based on the 
. principle of detecting radiation transmitted by a source of charged parti- 
cles and photons. Integral (10) has been calculated, for selected values of 
the parameters, by means of the everywhere convergent series as given by 
Hubbell et al. [9] or by more computationally tractable, but empirical ap- 
aaa the latter devised in particular for engineering applications 
[8]. 

Kalla et al. (10) have defined and studied the more general radiation 
integral 


2 


b 

a,b,c, A oa d = - . 

A ee | = 2 fi (x? +e) °F, («,6:n-35-) dz, (11) 
0 


where y > 8 > 0; a,b,c > 0;A > —1 (if b 4 oo, then —1 < A < 1) and 
2F;(a, 8; 7; x) is the Gauss hypergeometric function. Observe that 


H ew hehe: (12) 
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In Galue [5], Kalla, Galue, Kiryakova [12], Galue and Kiryakova [7] a 
generalization of the Hubbell integral and of (11) is studied in the following 
form: 


a,B,y 


b 
oa La 2\ B a? 
< An t A (a? +c) % (1-5) oF; (a.6:15-s) dz, 
0 


withe >0,0<a<b<o, wp >-l, A >-1, and if b > oo, then 
b 
-l<i\< 2a — 2-1. If p = 0, the function I reduces to H i oa 


a, 8,7 
defined by equation (11). 


In terms of the Appell’s double hypergeometric function, this integral 
takes the form 


ee = in TG TORE af | ares (ie ) 


x tP'(2? ++ at) “dt da . 


fan | 


(13) 


(14) 
From (14) and using the approximation formula 
(L+2/r)°! ~ (1—2)-(@-9/ (15) 


which in the neighbourhood of z = 0 is a better approximation to (1 + 
z/r)*—' than just the unity ({13]), in [6] we have obtained a single term 
approximation for I: 


Tel bee _ 74 pati TAH) E(y + o& +1) 
7 a, 6,7 ~ 8x c& r(4#8 + p+ a%) 


16 
Ny - caer aZI(y) 


rovided a?b*/(2c)<<1. 
Ty + aZ5)E(7 - 8) ae 


? 
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“In view of the symmetry of the hypergeometric function, equation (16) 
can be written alternatively as 
a,b,c, A, A) oa OH TAB Tu +a% +1) 
a, B, 7 _ 81 c& C A+3 b2 
eats a 
Tiy=a+ 6Z)0) 
Ty + BE)T(y - 2) 


4 


If in equation (17) 4 = 0, we obtain a single term approzimation for 
a, b,c, X 


namel 
a,By,y |’ : 


the function H | 


2 
pay oho soa PICS leet) 
Qa, B,Y 87 c& (ats +a%) 


Dy eho = 2) 


Then, the corresponding formula for the Hubbell integral is: 


oa = b T(8)r(e?+1) 
b)n~ — eZ 212 . 
h(a, dn Dea? TE+E) with a°b"/2<<1 


Following a similar procedure, we can obtain an approximation formula 
for a more general radiation integral. 
| a, b, c, 


e Single term approximation for the oe) a, (ap), (By) 
XP) \Y a 


function. 
Saigo and Srivastava [18] have studied a family of integrals of the form 


a,b,c, | 


_ @(P,9) 
a= Ree 
(18) 


ay 


b 

ap hed A pm —a ae gay py: |. 

_ An J: ve ma| Br, +) Ba; zm +¢ ae 
0 
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where min {a,b,c} > 0, € (—1,1) and pF, is defined by (1). 
One immediate special case of (18) obtained when p = g = 1 and 


m= 2, 1s 
(1) { @,,c,A] _ a,b,c, A 
hep bee | = ce ‘ (19) 


In [6] we have derived a single term approximation for SPP) =q)in 
the following form: 


b,c, 0 ga ptt 
S= (P,P) . age (ie a Ne m\—a 


r (ag) (2-2) wa (oop (He) Fa) 


xX ——— 


P (4H +08 +1) io P (Bisa — ap (25x) ) Pas) 


for 0<a™/c<<1(or a™/e>>1) and (ab)™/c? << 1. 
Letting p = 1 and m = 2in (20), we obtain for ap = a, a, = f, fi = 7, 


another single term approzimation for the function H . He ‘ in the 
form::. ie 
a, b,c, X ca bt! 2\—6 
re |= 8a cB cee)" 


2048) r(08 +3) r(o-0()) 1) 
r(eS+eE) P= 0 (ei) Me) 


for 0 < a?/c¢ <<1(ora?/e >>1) and (ab)?/c? << 1, and as a corollary, 
the corresponding formula for the Hubbell integral: 


von aet_6_ Heprte-rs) (0-4 (of) 
a, 8m (1 + a2)? I(2 + 82) re — 3 (s#)) 


for 0 <a? <<1 (ora? >>1) and (ab)? << 1. 
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The above single term approximations for the functions S and I are 
sufficiently accurate to allow a semi-quantitative assessment of the depen- 
dence of these functions on the parameters, which is usually an important 
first step in the physical applications. 


3. Single term approximations for the Gauss type 
generalized hypergeometric functions pi Fp 


Using successively the result [3] 


a*—"(1 — 2)?" Fy ((ap); (bg); we) de (23) 


i Oe 


B(a, 8) p+ilg+i ((ap), 0; (bg), @ + B; w) 


(a,Rea,Re® > 0;p < ¢g+1), Kalla and Galue [11], Kiryakova [14] ob- 
tained for pFy, p= q+ 1 the following multiple integral representation: 


F ij rere) = 1 1 1 
tee e By, »-5 Bp; . Ba, By —a) B(a1, Bz — a1) B(a2, Bs — a2) 


‘ 14-4 1 
<e ——— — hd l—w a1, A:-a-1 
B(ap-1, Bp — Op—1) LTA geo ) 
0 0 0 0 


x (1- wi JO aye 1 (4 — we)? hse so (1 — wp—1)%P717? 


x wer ora (1 —2(1 —w)(1 — w1)(1 — w2)-+ (1 - Wp-1)) 
x dwp—1 see dwy dw, dw. 
(24) 
The approximation formula 


1 — a(1 — w)(1 —w;)(1 — we) +++ (1 — wpa) & (1-2). + —) 
l-¢«z 


x (1+ —wi)(1 + wa) +++ (1+ -——wp-) for 0 < |z| << 1, 
—-2 _ 


(25) 
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and (15) turn representation (24) into 


ptil'p & Boe 1] = : Z : 


Bixee Bet B(a, fi — @) Suet B(a2, 83 — a2) 


111 
1 


a cery nrmey a | wen sane 


x wAr-8-1(1 _ 1)%1 +0 He z—)— 1 whe ay i jeteal ae Ly he—a-1 


X-+(1 ~ Wp) °P-2 F(R a)- Lyrae dwp-1:::dw2 dw; dw . 


(26) 
Then, from the known rational approximation of pFy (see Luke [177], 


85.12, (1)-(4),(8)), the following single term approximation can be ob- 
tained: 


: P11 (a. ie . 
p+ilp re fh ee 7 ~(l1—z)% I] z ( zs op( 72 5)) T(Bi41) 


1) wey Bp; i=0 r (Bix so ap(=*=)) P(a;) 


(27) 
with 0 < |z| << 1,Re (a; + ap(=4,)) > 0, Re(Aj41) > 0,7 = 0,1,2,..., p— 
1. 


Particular cases: 
a) For p=1 we get for the Gauss hypergeometric function: 
: + B(-))T 
- es = (1-2) F eee al eee (28) 
7 D(y + 8 (z25)) Pla) 


for 0 < |z| << 1,Re (a+ 6 (;25)) > 0, Re(y) > 0. 
b) lf p= 2, 


3, ee ~ (1-2) % 
D (6; + a (5) T(a) TP (62+ a2 (=) T(ai) 
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for 0 < |z| << 1,Re (a + ay (;2>)) > 0,Re(f,) > 0,Re (ay + Q2 (;25)) 
> 0, Re(62) > 0. 


4. Single term approximations for the Appell’s functions 


For the Appell’s hypergeometric functions of two variables Kalla and 
Galue [11] obtained the following single term approzimations. From (9), 


bc , By 
T(7) (a+ f+ ey) 


Fi(a, B, . ;2, ~ (1—2)7F(1—y)F 


(30) 


U 
with Re(y) > Re(a) > 0,Re (a+ ee fs) > 0. From (7) and (8), 
assuming 
(1-—az—by)~(1-—az)(1—by) if 0< |zy| <<], (31) 


and using (15), we find similarly, 


F(a, 8,8 57,73 2,y) & (1-2) (1 —y) 
Ely) Poy) £06488) T+) (32) 
PB) T(6) P+ 85) vy + 24)’ 


with Re(y) > Re(@) > 0, Re(7’) > Re(/’) > 0,0 < |ry| << 1 and small 
a; respectively, 


Fi(o, B,837,732,y) &(1—2)"%(1—y)~* 


Bg 
Pa) Por) Pb +5) r(a+é4) (33) 
FA) Ta) Ft) p(y fx) 
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with Re(y) > Re(8) > 0, Re(7 ) > Re(a) > 0,0 < |zy| << 1, small B. 
Analogous formulas hold for F3 and F4. 


Other single term approzimations can be established using (15) and 
(31) again, for example: 


' Tr I(y-a-fr-p 
Py (05.6) 8 39329) ACE CSS TT (34) 
provided. Re(7) > Re(a) > 0,2 ~0,y ~ 0, and 
F,(a, 8, B'37,732,Y) 
~— Ply) My - 8 - az) Ty’ - 6’ - ay) 
My - A)P(7')- 18’) Ty - ae) T'(7' — ay) (35) 


provided Re(7) > Re(@) > 0,Re(7’) > Re(6’) > 0,0 < |zy| << 1. For 
the technical details see Galue, Kalla, Leubner [6]. 


5. Approximations of generalized hypergeometric functions 
via generalized fractional calculus 


Formula (23) ({3]) turns to be also a corner-stone of the new and uni- 
fied approach to the generalized hypergeometric functions (1), Kiryakova 
(15],[16]. All the »Fq-functions are proved to be ”generalized fractional 
integrals or derivatives” of three basic elementary functions, depending on 
whether p < q,p=q or p=q+1, namely: 

COSg—pyi(t), ze expz, «%(1—- x), (36) 

These new differintegral expessions can be used successfully to reduce the 

problem of the pFq’s asymptotic behaviour to that of the well-known func- 
tions (36). 

The generalized fractional calculus developed by Kiryakova [16] is 


based on compositions of the so-called Erdélyi-Kober operators of frac- 
tional integration (E.-K. fractional integrals) (see [19]): 


: 1 — 6-1 7 4 
Ty Ha) = [ A  taot) da 8>0,7ER, (37) 
(8) 


-9]- 


generalizing the Riemann-Liouville fractional integrals of order 6 > 0: 
r ( je 1 r (1 ‘a 1 
eC ne pee seg PIED pete pe 
Re: fe). = / 7) f(t)dt =a i 7) f(ao®)do = x° I,” f(z). 
0 0 


Let m > 1 be integer, 8 > 0, %1,...,¥m and 6; > 0,...,5m > 0 be 
arbitrary real numbers. We define multiple (m-tuple) Erdélyi- Kober (E.- 
K.) operators of integration of multiorder 6 = (6),...,6m) by means of the 
integral operators 


1 
+ 6,)" 1 
1g g(a) = | Gini G ae | Hoot) ao, (38) 
0 


where the kernel-function Gin is a special case (m = p = q,n = 0) of 
the so-called Meijer’s G-function (see [3],[16],[17]): 


II F(% —-s) [[TQ—a; +s) 


api [ole gt |= [ei es 
re, AO Ae a Gesbeaeays TP ans) 
k=m+1 j=nti (39) 


Then, each operator of the form 


Rf(z) = wpb POE) PR) #15) with arbitrary 6) > 0, (40) 


is said to be a generalized (m-tuple) fractional integral. 


It turns out that single integral operators (38) are equivalent to the 
multiple integrals (like these appearing in Section 3, (24)) 


TO 7G) = ie aes eae iG 61 f(e)) } = i pe f(x) 


k=1 


1 Lem bp-1 Gk 
= (i= op) e~ sae 7 A Z a 
Sof | {I aa SC 4 20; 10h) do... dom 
o m 7 


0 


(41) 
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representing compositions of finite number m > 1 of commuting E.-K. op- 
erators ee k =1,...,m (Th. 1.2.10, [16]). This fact explains the name 


given above to (38). The Meijer’s G-functions are generalized hypergeo- 
metric functions of very general nature; they include as special cases al- 
most all the elementary and the special functions of mathematical physics. 
That is why, dealing with the single integral operators (38) instead of com- 
positions (41), has a number of advantages. The use of the simple proper- 
ties of special functions (39) (see [3],[17] 16, App.]) allows to develop a full 
chain of operational rules and theory of the generalized fractional calculus 
[16], based on the generalized fractional integrals (38),(39). The gener- 
alized fractional derivatives, corresponding to (38) are also introduced by 
means of explicit differintegral expressions, namely ((16], §1.5): 


Dein F(a) = Dy TEA” F(2), (42) 


where ny = 6x, if 6; is integer, or np = [6,] +1 if 6, is noninteger, k = 
1,...,m an 


aera in suitable functional spaces like the space of analytic functions, 
, 6 
Dee ( gee k f(z) = f(z). 


Both generalized fractional integrals (38),(39) and corresponding 
derivatives (42) will be called by the commmon name “generalized frac- 
tional differintegrals”. The following general proposition is derived in (16, 
Chapter 4]: 


PROPOSITION 5.0. All the generalized hypergeometric functions pFy 
can be considered as generalized (q-tuple) fractional differintegrals (88), 
(42) of one of the elementary functions (86), depending on whether p< 
@P=Gp=qrtl. 

To establish above proposition we use formula (23) stated in the fol- 
lowing form. 
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LEMMA . Let |x| < 00 (|z| <1 forp=q+1), then each pFy-function 
1s an Erdélyi-Kober fractional integral or derivative of a p—1Fg—1-function, 
namely: 


(ap )/1Cba)) pFq(a1, «++, Ap; 1, ..., bg; 2) 
Gy—1,bg—a : 
7 dae et Pp —1Fq—1 (G1, «+5 @p—13 b1,...,bg-1; 2)} af bg > ap, 
77 bg- —b : 
D,! me *"{p—1Fy-1(a1, vey Ap—-15 by, way bg—1; r)} af bg < ap. 


43) 


The three cases p < q,p = q,p = q + 1 are to be considered sepa- 
rately, first of them being more complicated and involving some auxiliary 
definitions and results. 


Ist case: p< q 


In {1] Delerue introduced a generalization of the Bessel function J,(z) 
for a multiindex v = (1,...,¥m), m2 1: 


(m) 3 (2/m+ Leek m, m41 
Jiy,.. Ym (2) = T(4 ie LL + 1) ofm ((v% a 1); ’ (x/m 3 1) ) ’ 
(44) 
called a hyper-Bessel function of order m. As a special case, the so-called 
(generalized) cosine function of order (m +1) follows: 


k a vate 1 k gh(mt1) 
COSm41(@) = oF'm (Ay) ;—(2/m+1)™ )= : en FT +1))!? 
(45) 


generalizing the elementary cosine function cosz =coso(z), m =1. 
In [2] Dimovski and Kiryakova proved a generalization of the Poisson 
integral representation of the Bessel function 


2 Gy spe? 1 
J)(«) = ioe) fo-e 2 cosat dt, v >— 5 (46) 


based on the Poisson-Dimovski transformation (see [16], Ch.3). In terms of 
the generalized fractional calculus, G- and pF -functions, the generalized 
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Poisson integral ([16], Th. 4.1.1, Cor. 4.1.4) can be written in the modified 
form: 


oFin (6, )3"5 -@) 


= ¢ LE I EM MH Cos gs ((m 4 Iet/im4y) 4 
== fom m[°| (aE) m y" a" cosms4i ((m + 1)(20)i/t49)) do, 
(47) 


where c = /(m+1)/(27)™ ES T(b;) and the condition b > k/ 
(m+ 1),k = 1,...,m is supposed. Otherwise, if some of the 6,’s does 
not satisfy it, its corresponding component in the generalized fractional 
differintegral should be considered as an E.-K. derivative ({16], Th.4.1.6). 
As an illustration, if b, := y%.+1 = k/(m+1)— 7, with integers n, > 
0,k =1,...,m, then (47) turns into a purely differential expression for the 
“spherical” hyper-Bessel functions (44) ([16],(4.1.42)), reducible in the 
case m = 1 to the spherical Bessel functions: 


(Ax 


1 
_ (22)42 dy {sos a 
—n-i ~ alae (dx?) - > V= 0, 1, Divs (48) 
Let us observe now that by p steps (23), 3 i.e.(43), a pFy-function, p < q 
can be reduced to a hyper-Bessel function, i.e. to a oFg— p-function: 


pig (Car) Ps (bi)t; ; a) 


Il Se I ace {oF (Cm ;2) \ » 


This ete differintegral relation combined with (47) (for m = 
q — Pp) gives the following particular form of the general proposition in the 
case p < q. 


THEOREM 5.1. The pFq-function, p < q 1s a generalized q-tuple frac- 
tional (differ integral of cosg—p41(z), namely: 
pFq(a1,.--, 4p; b1,..-, bg; —Z) 


= ATTA) { cosy—pis ((a —pt+ 12m) | , (50) 
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with A = v/q—p+1/(27)i-P i r(e)/TI rea and parameters 
j=l i=1 


k = 
na {E’ a{ PoE fork=1,...49—p 
Ak—qtp — 1, by — k-q4ps for k= q—pt+l,... (51 


If the conditions 


bp > ea ee 6 bo—prk > a, > 0, K=1,...,p (52) 


q—ptl’ 
are satisfied, then relation (50) gives Poisson type integral representations: 
pFq (a1,...,ap;b1,..., 69; -2 


1 
b,)? 
0 ( k k=1 

=A | oti > k 97? (a | 

0 Q—Pti kay? \WR—9+P/ kag—p4ti 
x 07 cosy—par |(q— p+ 1)(20) 7-9] a5 

rv fa (1—t,)* pA 

SL @ k 

0 0 = 


k q bd = 
. el Il (1-t,) k~Tk-qtp a 
(53) 


1 
X COSg—p+1 l(a —pt1)(xt,... tg) PF | dt, ...tg. 


If (52) are not true at least for one of the indices k, then y(%e)s(6e) in (50) 
is considered as a generalized fractional derivative of form (42). 


Relation (50) including integral representations (47),(53) and their dif- 
ferintegral analogues, expresses the pFy-functions, p < q as fractional dif- 
ferintegrals (integrals or derivatives) of the generalized cosine function. 
This fact generalizes the well-known representations like (46),(48) of the 
Bessel functions via the cosine and thus, suggests the name Bessel type 
g-h.f-s for the functions pFy with p < q. The Bessel function itself is the 
simplest special function of this class of g.h-f-s. 
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2nd case: p= q 
Applying relation (43) to a function pFp consequently (p — 1) times, 
we reach to a ;/4-function which on its side is representable as an Erdélyi- 
Kober operator of the elementary function 2% exp z: 


[ee-aret) ifb>a 


pret) ifb<a. 
(54) 
Thus, combining (43) and (54), by p steps we obtain the form of the 
general proposition in this case. 


I'(a) 


1—a 70,b-ay a-1 ,& 
T(b) 1F\(a;b;2)= 27 "T,” “{x »={ 


THEOREM 5.2. If p= q, the g.hf. pFp(x) is an p-tuple fractional 
integral or derivative of the elementary function {r“—1e"}, namely: 


Pen ras De naibes = Pi gia {Teed {z%—1e%} (55) 


Pp 
where yp = ap—a1, dp = bp —az, K=1,...,p and lr’ = TT (P'(3;)/T(a;)I. 
j=! 


If 


b. > a, > 0, & =1,...,p, (56) 
this relation yields the following integral representations: 


i 
b 
pFp(a1,..-, ap; b),...,b9;2) =T | o88 0 ey Jom ' exp(xo)do 
0 


1 — tp)Pk- 9K ee 1 
=i To fT a EG aap exp(at ...tp)dty ... dtp. 
57) 


The above theorem justifies separating of all the g.h.fs pFp ina as of 
g.h.f-s of confluent type, involving the confluent hypergeometric function 
1F;(a; 6; x) = G(a; b; x) as a simplest case. 

For parameters not satisfying (56), relation (55) gives differintegral 
expressions. For example, we can introduce “spherical” g.h.f-s of confluent 
type (by analogy with (48)) representable by pure differential operators of 
exp z. 
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COROLLARY . Let all the differences aj, — b, = np, k = 1,...,p be 
nonnegative integers. Then, the “differintegral operator” in (55) turns 
into a differential operator Dy of integer order n = m +--+ + np = O, 
namely: 


Pe (by) 
pp (b1 + m1,--- 5 bp + np; bi, ..., bp} 2) as iétap 


ja b; T nj) 
(58) 
x Il Hex +b, +7 -—1)| {exp z} = Qp(z) {exp z}. 
k=1j=1 


Differential representation (58) gives an example how differential for- 
mulas for the “spherical” g.h.f-s introduced in [16] can be used for their 
explicit calculation, especially in the case p = q in the form Q p(x) {exp z}, 
where Qp(zx) is a p-degree polynomial. 


8rd case: p=qt+l1 


The generalized hypergeometric functions pFy, p = q + 1 are said to 
be g.h.f-s of Gauss type (see [16], Ch.4) and are considered for |z| < 1. 
In this case the starting specific result is representation (6) of the Gauss 
hypergeometric function, written as an E.-K. fractional integral: 


I(a;) ‘ ’ 
T(bi) 2F} (a1, a2; bi; x) (59) 


= taf (a = gy) = gl fo 1(1 — 2) ™} 


for 6; > a2 > 0, or as a corresponding fractional derivative if az > b, > 0. 


Since by (q—1) steps (43) a p»Fg-function reduces to a 2F-function 
and the composition of fractional differintegrals in (43) and (59) gives a 
g-tuple integral or derivative, we obtain the third form of the above stated 
proposition. 


THEOREM 5.3. In the unit disk |x| < 1 the g.h.f-s of Gauss type 
pia, P= a+1 are q-tuple generalized fractional differintegrals of ele- 
mentary functions of the form #%(1 — ay, namely: 
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qtilg (Gryueiy gai? bya Dg) 
= r’ ia globe Da (Pam aeas {2™-1(1 = x)~% } ; (60) 
with TY = T]4_,(0(;)/T(aj41)]- This means that for 
by > ap >0, k=1,..,m (61) 


the following Poisson type integral representations hold: 


(by) 


1 
—p" q,0 
qtify (+e) =T [ot [> Ge 


| o(1 4 20)-"do 


0 j=l 


q 1 1 q : 

T(b;) iE {TI bp—aps, —14%k41—-1 

= 1 — t,)?h7 ok+1 yy ka 

LG rcreyccrerry 1 Caj4a PCG; ia 541) (q) II | ( k) k | 
(1 ates tg)” ies Loa 


(62) 


COROLLARY . For g = 1 representation (62) coincides with (59), writ- 
ten as the known Euler formula (6) for the Gauss function: 


PG (ft) erat ths —Ga—1yaa—1 
FE > by; 6 
aE ON ae) NT aa) | 7 reapag, o *AG) 


valid in |x| <1. 

This formula proposes a way for an analytical continuation of 2Fy(x) 
outside the unit disk to the domain |arg(1 — x)| < 7, where the right- 
hand side of (63) represents an analytical function of x. For the same 
reasons, formulas (62) can serve as analytical extensions of the g.h.f-s 
gtifg(£),g > 1 outside |x| <1. 

The case with parameters not satisfying condition (61) yields general- 
ized fractional derivatives in (60) and provides also useful corollaries. By 
analogy with the previous two cases, we introduce the notion spherical 
g-h.f-s of Gauss type when all the differences a, — by = np, k =1,...,q are 
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nonnegative integers and 941 F(z) is representable by a purely differential 
operator of a function (1 — z)?. 


Another interesting case concerns the so-called hypergeometric poly- 
nomials 


nr 

(—n), (a1), ++ (ap), ak 
Fg (—1n,ay,-..,4p3b1,...,bg3 2 =) ee 

purky ( , a 4) c+ (bigs (by), 


when p = q. By taking ag41 = —n, n > O-integer and a, > by > 0,k = 
1,...,q, the fractional derivative form of (60) turns into the Rodrigues type 
formula 


(64) 


T'(a;) 
ne pti ly (—n, a1,...,493b1,.-+,bg32) 


= per 1),(a,—by) {(1—2)"} (65) 
gpl a) {211 2 x)"} — pimbg ptq—bg tp bg—1 


x D%-1—9-1 a pts 02 Pyt2—bo phi pub ema _ x)" } . 


Some special cases of (65) yield classical Rodrigues formulas. For ex- 
ample: p = q = 1 with aj = n+1,b; = 1 ands —- ie yields the 
Rodrigues formula for the Legendre porgnoniars 

Pa(z) = (-1)"2Fi (mn +1;1; 5 x 
_(-1)" a fl-2"1l4+2"] 1 
nt dzt[ 2 2 
and p = q = 2 with qy =n+1,b) = 1,a2 = (62 = p(¢ >p> 0) gives the 
Rodrigues formula for the Rice polynomials, viz. 
Rn(z) = 3Fo(-—n,n+ : aC else a) 
a, T(p) 1—p C=p n (67) 
~ niT(C) ce Oe =) ue eee) 


All the results mentioned in this section show that there are, essen- 
tially, three kinds of g.h.f-s »Fy, similar in properties and reducible to the 


d” e (66) 
hana 2 ; 
~ Onn! dan (2 = 1) bi 
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three elementary functions (36) and to the three “initial” (simplest) g.h.f£ 
Ss of), 171, 2f;. Then, many results for them, including approximation 
formulas, can be transferred from the simplest functions to the general- 
ized hypergeometric functions, using that in general, the Erdélyi-Kober 
fractional differintegrals and their compositions preserve the asymptotic 
behaviour near the origin x = 0. 
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1 Introduction 


We denote the complex plane by C. If D C C is a complex domain then a 
rational map of the form 


az+b 


d 
f(z) rarer | ad—be #0, z=—-~ €D, 
is called a Mobius transformation. Using a matriceal representation of f, we 
shall study some properties of the iterates of f provided a,b,c,d € C are in- 
- dependent random variables. Applications to branching processes are given. 
Throughout this paper we will denote by F the set of all Mobius transforma- 


tions that is 


az+b 


di ere | 


d 
ad — be £0, z#——}. 
Our paper is organized as follows: Section 2 deals with the iterates of f, in 
Section 3 a generalization is given, Section 4 shows how it is possible to obtain 


fractals structures using Mobius transformations and Section presents numerical 
results. 


2 The Iterates of f 


a b 
Denote by A = ( ot 


of z + (az + b)/(cz +d). Consider e; = ( ; ) a= ( : ) , {e1, €2} being 


) a 2 x 2—matrix whose elements are the coefficients 


the canonic basis of C and Z = ( : ) , be acolumn vector. We remark that if 


f © F then 
f(z)= (ce AZ)/(e5.AZ) (2.1) 


where we denoted by z* the transpose of z. If on F we consider the composi- 
tion then F becomes an abelian group. Moreover, if f,g € F have A and B 
respectively as associated matrices then 


(f 09)(2) = (e, ABZ)/(e} ABZ) 


hence (F, 0) and (Moxa, :) are isomorphic. This fact allows to write the iterates 
of a function f € F as follows 


f(z) =(f ofo...0f)(z) = (eh A"Z)/(e,AZ). 


In this manner, the study of f() is transferred to the study of A”. 
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If we denote by Z, = ( o ) then the iteration formula may be written as 


follows 
Zn41 = (e} A" Z1)/(e$A"Z1), ne N* 


We shall use Perron’s formula applied for the n—th power of a r xX r—matrix 


k F 
1 d™i-i 
A? = S> ——___ a MOS WMO aAy 
= (m; = 1)! d\mj-1 { ( i) ( ) } ee 


where m1,..., mz are the algebraic multiplicities of eigenvalues A1,..., Ax, ™1+ 
1. tm, = 7, k <r. In our case, if Ay, A2 are the zeros of the characterictic 
polynomial of A, A? — (a+ d)A + ad — be = 0 then 


poo t+ tee (+ b i) 


If A1 = piexp(ia1), Az = prexp(ta), pi,p2 > 0, 21,2 € [0,2m) are the 


complex eigenvalues of A, then the iterate f() may be written as 


72) = [ott Eas — pot Ena — 4 (oR En — pp En)] 2 +b (op En — ph En) 
e (pRB, — py bz) z+ pit Ens, — p27 E24 — a (pf Ep — pp ER) 
where we denoted E* = exp(iagn), k = 1,2. 
Taking into account these general backgrounds, we shall study the following 
cases: 


n 
I. pi # p2, % F a2. Suppose pi > po and denote r, = (2) . After some 
straightforward calculation we get 


lim f(z) = [p1 exp(ia1) = dj z+b 
n-+00 ez + pi exp(ia,) —a - 


But A; is eigenvalue of A, therefore 


lim f(z) = ae é 


A - a 


Hence, independent of zo the sequence of iterates { f™(20)}, 
to the complex number (A; — d)/c, c# 0. 

II. a,b, c,d € R. If we denote 1 = pexp(iw), A2 = pexp(—ia), with p > 0 
and a € (0,27), then 


>9 18 convergent 


[psin(n + 1)a — dsinnalz + bsinna 
(esin na)z + psin(n + l)o — asinna 


f(z) = 
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Suppose that exists k € N* such that a = §. Ifn = kp, p € N* then 
f(z) = z, that is the identity of F. For n = kp —1 we have f((z) = 
(—dz + b)/(cz — a) = f-1(z). Therefore, the iterates { f()(z)} have the form 


f(z), FA(z),..., FEM) =f), fHO(2)=2 
fFO+D(z),  fE42)(z),...,  fAR-D(z) = fol(z), fCH(z) =z,... 


We can deduce that the sequence z, = f(z), n = 0,1,2,... has a finite 
number of points (cyclic points). These points belong to the circle of equation 


2% + ay(z + Z) — iao(z — 7) +03 = 0, 


and @1, a@2,a3 € F are solutions of the following system 


24 — 21 i(z _ 21) 1 ay —|zi|? 
29 —~ Zo i(Zo - 29) 1 a2 = —|z2|? (2.2) 
23 — 23 i(Z _ 23) 1 ag —|z3|? 


where z; = f(z), 22 = f(a), 23 = f(z2) with zo fixed. If does not exist k 
such that a = , then the iterates f()(z) belong to a circle of equation (2.2) 
and the sequence {f™(z)} 5, with zp fixed is mapped into the points of this 
circle, [3]. 7 

III. a,b,c,d € R,.’ According the theorem of Perron-Frobenius, a positive 
matrix has a dominant positive eigenvalue and the corresponding eigenvector 
has positive elements. Let’s denote this eigenvalue by A; and let us consider 
A2 < A1, A1, A2 € R. We denote by ui (v1), ue (v2) the corresponding left (right) 
eigenvectors. From the normality condition < u,;,v; >= 1, i= 1,2 where <,> 
is the inner product we have 


A” = Mozui + AZ vgus. 
In this case, after some calculations, we get 


(At m€1 + AZ 01 P1) z+ AP me2 + AF Po 
tmes + AZ a2/1) z+ AT 1 2€2 + Ab 280 ; 


Here we denoted us =(7 ),u=(%),u=(%), m=( 5 ), 
2 2 2 


and m; > 0, ¢; > 0, 2= 1,2 satisfy m1, + m2€2 = 1. Dividing by A? and passing 
to limit we have 


f(z) = 


m f(z) = 4 
jim f (2) 72 


hence in this case, independent of zo, the iterates are convergent to a point 
located on the real axis. 
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Application. Let’s consider a branching process, [1], {,,n € N} whose gen- 
erating function is given by 


l—b-e bz 
as ere eee 


0<b,c<1, b+ce< 1. 


f(z) is the sum of power series )°,5, pez”, with py = be*-!, k = 1,2,... and 


Do = ates, In a branching process, £, represents the number of individuals 


from a given population at n—th generation and the expectation will be 


m= f'(l)= >> kpe = 


nye 
Post (1 —c) 


The matrix of f will be 


2 = 
gat c 1-b .) 


c“-c¢ l-ce 


and its eigenvalues are 4; = b, Az = (1—¢)?. 
From (2.4) we obtain 


lim f(z) = 


n—-oo 


Am +6-1 
e—1 ’ 


where Ay = max{A1, Az}. If m > 1 then f(z) > 1='5¢ while if m < 1 then 


f(")(z) 1. The limit of f()(z) represents the probability of extinction of the 
process, [1]. 


3 Generalization 


Let A = {Ai,A2,...,Ar} be a set of r 2 x 2— non- singular matrices which 


etre . b 
correspond to some Mobius transformations form F. If A, = = . ) E 
n 7 


A, n> 1 then we build the sequence as follows: 


z given, 21=fi(z), 22 =(fe° fi)(20),..-- 


Taking into account (2.1) we obtain the matriceal representation of this sequence 


as 
ef (An...A1)Z Zo 
= Z= 
PE Ay elie: ee 

Let us consider W = C, X = A, W the borelians of W, and 4 the set of parts of 
X. Consider a probability measure {p,(w), w € Wh, cx with pe(w) = pz > 0 
and )>,¢x Px(w) = 1. Moreover, for all w € W we consider the family of 
applications f:W x X — W which are (W @ 4, W)—measurable, f € F. The 
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tuple {(W, W), (X, ¥), f, p} is called random system with complete connections, 
[4]. 

We remark that the sequence obtained above, (zn),5 9, where for each n a 
random f,,(z) is choosen, defines a Markov chain. The transition probability of 


this Markov chain is 
P(z,A) = > La(fo(2)) Pe (3.1) 
TEX 


for all z € W and A € W. This is the Markov chain attached to the random 
system with complete connections. 

Let us suppose without loss of generality that the matrices A, are posi- 
tive (the general case is similar). For practical reasons we have to study the 
asymptotic behaviour of the following random product of matrices ° 


A) = A,...AgA1. 


For this purpose, we consider A = ( : ; ) a positive matrix and its decom- 


position as a product of a diagonal matrix and a stochastic one as follows: 


2preh 24 25( gaa eer 
o=("S) ota) PoC Be). 


Therefore A = D- P. Using this decomposition for A‘) we have succesively 


Al?) = Dy PnAn-1 id AY = Di AnatAans on AL 
= Dy Dn-1Pn-2An—-3.--A1 =+::=Dy,...DnPn. 


Here Di,...,Dn are diagonal matrices, P, is a stochastic one and An_1 
P,An—1. In general, if we denote by A the product between a general matrix 
and astochastic one that is A = AP from the componentwise interpretation we 
get 
mina;; < mina; < maxa;; < maxa,;. (3.2) 
aj a,j a,j a,j 


If we denote 


La) or) nm) 9 
vs ( fm) qn) , D-Day = ( 0g) ) 
- Pin 1— pin 
a ( 1— pon P2n ) 


A) = ( AO pin AL”) (1 = pin) ) . 
g™) (1—pon) gan 


Based on the inequalities (3.2) and taking into account the positivity of Ap, it 
results that exists 6 > 0 and ng > 1 such that pin > 6, pon > 6 for alln > no. 


then we have 
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3.1 Proposition. In the conditions described above, there exists and it is finite 
the limit 


lim id log a”), 
n—-00 TN 


Sketch of proof. Indeed, we have 


A Joga = Llogh™pi, = Llogh™ + + log pin 
= LVL, log hi + 4 log pin 


where A(") = hy---hy with hy > 0 for every i = 1,2,... But for n large enough 
Pin > 0 and then + log pin — 0. From the law of large numbers it results that 
lim 4 )77_, log hj exists and it is finite. Let’s denote this limit by h. Therefore, 
+ log a‘) —+ h. For the other elements of A") the demonstration is similar. 

By Proposition 3.1 it results that if n is large Saye then A‘) approaches 


e” 


ong eng , where h and g are 


to a matrix with proportional rows A“ = 


positive random variables. Therefore, the sequence (zn),59 is given by the 
recurrence A ‘ ~ 
n m™m 
Zz ™ € td + € — n(h—-g) 
n+ e 
eroz + end 


hence independent of zo. 


4 Fractals from Mobius Transformations 


Let (C, D) be the complex metric space endowed with the distance d(z1, z2) = 
|z1 — z2|, 21,22 € C and {f;,..., f,} Mobius transformations from F whose 
matrices are A;,..., A, respectively. Let us consider also the set 


» Po >; YD) pe =1. 


TEX 


{pe,2EX} X={I1,2,...,r} 


It is obviuos that for all k € {1,2,...,r}, f, are Lipschitzian complex functions 
defined on C, that is there exist s,; such that for all z, w € C the inequalities hold 
d(fx(z), fe(w)) < sed(z, w). Barnsley, [2], proves that the Markov chain (zn),59 
built in Section 3 having the transition probabilities (3.1) and the transition 
operator 


(Ta) (2) = f a(y)PCe,d) = Sealeo 


has a unique invariant measure without the eile that f, would be con- 
tractions and the space of states would be compact. 


1.0 


se 


Genenenesenr nests sneees fii mare mam 
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4.1 Theorem. [2] Suppose that there exists 0 < r <1 such that 


d(fi(w); fe(z)) |" < 


d(w, z) (4.1) 


uniform with respect to w,z €Y,Y being a complete metric space. Then, the 
Markov chain (zn), 59, built as above, admits a unique invariant measure }. 


The support of pw is a fractal. In our case, the states space is C' and for a 
set of Mobius transformations f),...,f, there exist pi,..., Pr; py per = 1, 
pr >-0 such that (4.1) holds. 


5 Numerical Examples 


First, we remark that any Mobius transformation may be written as 


joys 


Zz+Y7 

It generates a random number q. If q € (0,0.8) then a = —3+ 5i, 8 = 0.8 —3, 
and 7 = —4+ 51. If not, for g € [0.8,0.86) choose a = 2+ 0.41, 8 = -0.7-i, 
y¥ =1-2ielse a = —44+4i, 8B =0.7-i, y = 5427. Compute and plot arandom 
orbit of zo under f. In Fig. 1 is displayed the result after 150,000 iterations. 

A special case of iterative relation with applications to n-dimensional branch- 
ing processes is 2441 = f(zz), k > 0, cy € R”. Here f(x) = (2,... 2K?) 
has the components of the form 


5 6 P, VEC. 


ajay +a?ag+---+ aan +a? 
byay + bata + +++ + dpa” + bo’ 


Ulam’s conjecture, [5], states that the points z; belong to a cone in R”. For 
n = 2, we have the following iterative relation 


fi(z) = 


7=1,2,....n, @= (@1,..., 2). 


eee ai, + aby, + a} eles azz, + azy, + a2 
fend by xp + boy, + bo’ + biz, + boyy, + bo 


If we choose two iterative relations of the form (5.1) which succeed each other 
with probablilities p and 1 — p, 0 < p < 1 respectively then we find a Markov 
chain. For 150,000 iterations, the suppurt of the invariant measure of this 
Markov chain has the appearance given in Fig. 2. Here we have: with 0 < p< 
0.5 


(5.1) 


nes dr, + 3ye —1~ eee tet lye tl] 
ate Fee oes OS Die, hay eS 
else 
1.50, + 2.ly, + 1.2 2c, + 1.6y, — 2 
Thi = Yeti = 5 


“2.4, +2.8y, +2” ~ 2 Ar, +2.8y, +27 


So, Ulam’s conjecture is confirmed. 


-112- 


References 


[1] K.B. Athreya, P.E. Ney, Branching Processes, Springer-Verlag, 1972. 


[2] M. Barnsley and J. Elton, A New Class of Markov Processes for Image 
Encoding, Adv. Appl. Prob., 20 (1988), 14-32. 


[3] P. Halmos, Lecture on Ergodic Theory, Chelsea Publishing, 1956. 


[4] M. Iosifescu, S. Grigorescu, Dependence with Complete Connections and 
Its Applications, Cambridge Univ. Press, 1991 


[5] S. Ulam, Problems in Modern Mathematics, Interscience Publisher, N.Y., 
1960. 


ALGEBRAS, GROUPS AND GEOMETRIES 40 113 — 133 (2024) 
DOI 10.29083/AGG.40.01.2024/SC113 


A SERIES-ITERATION METHOD IN THE THEORY 
OF ORDINARY DIFFERENTIAL EQUATIONS 


D. Dimitrovski 
Institute of Mathematics 
Faculty of Science, Skopje 
M. Mijatovic 
Institute of Physics, Faculty of Science, Skopje and 
Institute of Nuclear Sciences 
Vinca, Belgrade, Yugoslavia 


Abstract 


The main contribution of this paper is the treatmnent of the the- 
ory of linear ordinary differential equations by the well-known traditional 
methods of series and iterations. In the iteration method, instead of the 
continuity we require the analyticity of the coefficients of the differential 
equations. Then the iteratioris can be carried out for the arbitrary form 
of the coefficients. Our procedure determines the fundamental system 
of particular integrals {y1,y2,...,yn} with arbitrary precision. There- 
fore our theory of differential equations is non-Liouville. (In the Liouville 
theory we know only that the solution has the form y(z) = S*p_» Cee, 
without a procedure for finding the particular integrals). The Wronskian 
W(y1,Y25+-+) Yn) = 1!2!-+-(n — 1)! of the fundamental system of a differ- 
ential equation of n-th order has a new meaning. It becomes the basic 
metric relation in the space of particular solutions. Knowing that linear 
differential equations of the second order with constant coefficients corre- 
spond to usual trigonometry there exists a natural generalization of the 
trigonometry (elliptical and hyperbolical). From the Wronskian as a basic 
trigonometrical relation the Euler formula, addition formulae etc. follow. 
The solutions in the form of a series of multiple integrals provide a nat- 
ural approach to approximate calculations. The approximate solutions 
are especially suitable in the case of monotone coefficients of differential 
equations. Working in the space of the fundamental particular integrals 


¥1 @ y2 ® +++ @ yn one can obtain interesting geometrical relations. 
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1. The Basic Principle and Definitions 


The expression 


y™ 4 ay (x)y"-)) + ap (x)y("-?) +++++dn(z)y = b(2), (1.1) 


where a1(z), a2(x), +--+, @n(x), b(x) are analytical coefficients is called 
linear ordinary differential equation of n-th order. 

In the theory of differential equations the following well known prin- 
ciple is valid: 

The solution of a differential equation depends of the coeffi- 
cients on the equation only. 

This principle cannot be proved because it is not possible, for a lot 
of differential equations, to express the solution via the coefficients. Some 
equations are not solvable, some have multiple solutions, while of some is 
known that these have solutions but the algorithm for their construction 
do not exist (only the theorem of existence). However, practical experience 
during the history of differential equations shows that the above principle 
is valid. 

This paper is excellent ’proof” of this principle. We will to try to 
show that all types of solution (general, particular ...) depend uniquely 
on the coefficients of the equation. 

Definition. The solution of differential equation (1.1) is called a 
general quadrature if it has the form of a series which can be expressed 
by the coefficients of the equation (1.1). 

We now enumerate three different possibilities: 

a) The solution has a closed form in terms of well-known functions, 

b) The solution has the form of a series, not of z but of the coefficients 
of the equation (1.1), 

c) The solution is a power series where the coefficients depend on the 
coefficients (1.1) explicitly. 

We will call all these possibilities quadrature solutions in the wider 
sense. 
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2. Linear Equations of I Order 


First we consider the simplest possible case, the linear canonical dif- 
ferential equation of the first order 


y' +a(x)y = 0, (2.1) 


with the intention of explaining the series-iteration method. 
The first step in our consideration is to integrate differential equation 
(2.1) and get 


y = —fa(a)y de +C, (2.2) 
0 


where C' is the constant of integration. The next step is to substitute the 
all right hand side of the equation (2.2) instead of the variable y in the 
right hand sade of the equation (2.2) 


= ere [- fotery dx+C| dx +C, (2.3) 
0 


and we obtain 


y= fala) de t a(x) y dz — Cf a(x) de +C. (2.4) 
0 0 0 


The next step is to substitute the variable y from (2.2) in the right 
hand sade of the equation (2.4). If we continue this process of iteration 
we obtain the following form of solution of the differential equation (2.1) 


y=C |1—fa(2)de + f a(x) de f a(x) de 
ae eres (2.5) 
= Jala) de f ale) de f a(x) de + | | 


0 


Using the identity (Dyson formula [1]) 


zr Zz Zz zx k 
f a(x) de fala) de fale) de = all a(z) ae) (2.6) 
ey 

k 
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the solution (2.5) can be transformed into 


y=) a (- jo(e) de 7 (2.7) 


k=0 


which is, of course, the series of the exponential function, i.e. 
y = Cem fo o(2) de, (2.8) 
Similarly, the nonhomogeneous differential equation of the first order 
y! + a(z)y = (2), (2.9) 
has a first integral 
zx x 
y=—Jfa(zr)ydze+ fb(x)dr+C. (2.10) 
0 0 


Again if we substitute the all right hand side of the equation (2.10) instead 
of the variable y in the right hand sade of the equation (2.10) we find 


y=—fa(z)|—fa(x)yde+ [b(e)+C| de+f2(a)de+C 
0 0 0 0 


= f a(z) dx i a(x)y dx — f a(x) dx [ b(2) dx — Cf a(x) dx + i b(x) dx + C. 
0 0 0 0 0 0 


(2.11) 
Continuing step by step the process of iteration we obtain the solution 
of nonhomogeneous differential equation (2.9) in the form 
yo Ce tome) dey. (2.12) 
where Y, is the particular integral 
i So So oF x So 
Y, = J b(a) dx — J a(x) dz f a(x) dz + J a(x) dz [ a(x) dz [ W(x) dx —--- 
0 0 0 0 0 0 
(2.138) 
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After differentiation of the expression (2.13) we have 


VY) =) 0s) f 6(x) de zs [ a(a) de J (2) ds decealic, - (0aa) 


So, the particular integral (2.13) satisfies the identity 


Y,' + a(z)¥p = b(2). (2.15) 


Because the identity (2.15) is the same as a differential equation (2.9), 
the particular sulution Y, is equal to well-known particular solution of the 
differential equation (2.9) 


rae fo o(2) dx [d(x)e! a(x) da oa (2.16) 
0 


The main attraction of the method of iteration is that in the formula of 
solution (2.12) appears naturaly the solution of homogeneous part (with 
a constant of integration) an also the particular integral. 


3. Linear Equations of I] Order 


Let us study the canonical differential equation of the second order 
with one coefficient [2-5] 


y" +a(x)y =0, (3.1) 


where a(x) is an analytical positive function in the interval [z1, 22]. Equa- 
tion (3.1) is well known and with a special form of a(z) it becomes the 
equation of Mathieu, Lamé, Hill, Bessel etc. In classical mechanics equa- 
tion (3.1) represents harmonic oscillator with time-dependence mass. Also 
this type of equations is equivalent to stationary Schrédinger equation 
which plays a fundamental role in quantum mechanics. 

Now we applied the method of iterations on the differential equation 
(3.1). After the integration we have 
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y=— fa(e)y dz + C2, 
Sz (3.2) 

y=-—fdzfa(x)ydex+Cor+Ci, 
0 60 


where C} and C2 are the constants of integration. Repeating the procedure 
in the previous paragraph we substitute the all right hand side of the last 
equation instead the variable y in the right hand side of the same equation 
and we find 


a(x) da? —.. | se 


fax) de? [va(x) da? =e | : 
0 0 


So, the solution of the differntial equation (3.1) has the form of summa 
of two series of double integrals 


y(z) = Ciyi + Cryo, (3.6) 


where 


< 
ee 
Il 
M 
| 
fry 
wa 
> 
ome 
ome 


a(x) dx? ff ale) de Pf Fale) de®, (3.7) 


—— 
k double integrals 
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=a2+ S\(- i J fa(e) ) da? [fala a) dz Pf eae) dx? . (3.8) 


k=1 REST” 
k double integrals 


The form of solution as a summa of multiple integrals is good possi- 
bility to use some software for analytical integration. In the reference [6] 
we have used the Mathematica [7] to perform integrations in the case of 
a(z) = {e*,]ln2z,cosz,sinz}. 

Theorem. Let a(z) be a analytic coefficient (not necessarily pos- 
itive) in the interval (0,X], then the series (3.7) and (3.8) converge uni- 
formly for all z € (0, X]. 


Proof. Ifa(z) is an analytic function it is continuous and bounded. 
So there exists the number M such that |a(z)| < M4. Then we have the 
majoranta 


<14MJ J a? +M? ff da? f f de® +--+ (3.9) 
0 00 
x? sa a4 yen 
=1 —+Me-—+4+.---+ MM” es 
tM AMG ti My + 


oe 2n 
nl< ce (3.10) 


n=0 


Using the D’Alamberts ratio test, for the n-th term in the series (3.10) is 
valid 


F,.\2(n +1) 2 


is atl oy (IE =p es 
oe An = (Mz) =e n+1 ai ist) 
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Then the series (3.10) converges for any z. Therefore, the smaller se- 
ries (3.7) with positive and negative terms converges uniformly for,any 
analytical coefficient a(z ) QED. 

Similarly, for the series (3.8) we can write 


|x| dx? +... (3.12) 


= |2 saat ? i +- 


OF vy 
(3.13) 


bal < OM a 2 


Because the last majoring series, according to the D’Alamberts test, con- 
verges, the series (3.8) converges uniformly. 

- The consequence of the previous theorem is that in the series (3.7) and 
(3.8) the usual analytical operation of such differentiation and integration 
term-by-term can be applied. 

Let us consider the homogeneous differentail equation 


y" +a(a)y' + b(z)y =0, (3.14) 
where a(z) and b(2) are analytical functions. With substitution 
= en 2 fo s)dz, , (3.15) 


the equation (3.14) can be transformed into a canonical type 
z" 4 A(z)z =0, (3.16) 

where 

a'(x) 


2 


2 


AG Sia sa(2) 2 (3.17) 


Let us we study the most general case of the differential equation of 
the II order, so called the nonhomogeneous equation 
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y" +a(z)y’ + W(z)y = f(z), (3.18) 


r), b(z) and f(z) an analytical coefficients. Then by the substi- 
15) the above differential equation transforms into 


where a(z 
tution (3. 


2! + A(o)e(2) = file), (3.19) 
where A(z) is determined by the formula (3.17) while 
faa) = 8 8% F(a), (3.20) 


So, instead the equation (3.18) we can to study the equivalent equation 
(3.19). 
Theorem. The differential equation (3.19) has a particular solution 


(3.21) 


= 
8 

ee 
QQ. 


rz uz 
2° TT (Ff Aleyaet) ate) de? - 
00 \00 
The proof follows directly by differentiation of the expression (3.21) 
and substitution into the equation (3.19). 


Theorem. The general solution (in wider sence) of nonhomoge- 
neous differential equation of the II order (3.14) has the following form 


> 
8 
p, ae 
= 


2? ff A(e ) dx? ane mai int —} 
00 
(3.22) 
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4. Generalized Trigonometry of Canonical Equations 


In differential equation (3.1) by putting a(v) = 1 we obtain the equa- 
tion of a harmonic oscillator. Then the solution (3.7) reduces to 


Ys 2 ctr rar 9 
yy =l—ffl-de*+ff1-da°*ff1-da°—.-- 
00 00 00 


4.1 
= a? gt -98 Gy) 
=o a. rien nay. 

and it is standard cosine function 
Y1 = COST = COSg(z)=1 2. (4.2) 
And analogously 
yg =uv—ffl-ade*+ff1-de*ff1-adz* —-.-- 
00 00 00 
ae a oe (4.3) 
Ger gt 
which denotes sine in usual trigonometry 
Yo = sin? = sing(z)=1 2. (4.4) 


Generalizing the above facts [3] we make the following 
Definition. The generalized cosine is the solution of the differential 
equation (3.1) (a(2) > 0) with the initial conditions 


y(0) =1 and y'(0) =0. (4.5) 


Therefore 
yi(@) = coSq(n4)® =1—f fa(x)dx? + f fa(x)d2? f fa(x)d2” —--- (4.6) 
00 00 00 


Definition. The generalized sine is the solution of the differential 
equation (3.1) (a(z) > 0) with the initial conditions 


y(0)=0 and = y'/(0)=1. (4.7) 
Le. 


Yyo(2@) S sing) 2 = 2— [ f we(x)dz? + / i) a(a)de? ff re(a)dz? —--+ (4.8) 
00 00 00 


Now we will study the derivative of the solutions cosg(z) # and sing(,) 2. 


Theorem. The generalized cosine (4.6) satisfies the identity 


Oo 
(cosa(.2) x)’ = — fa(x) cosg(,) 2 da. (4.9) 
0 


Proof. According to (4.6) we have 


(cosa(x) x)’ = —fa(x)dx + fa(x)dzx f f a(x)dx? — 
0 0 00 
=— fe(z) 1- ff a(a)dex? + [ fe(z)dx? i f a(x) dx? — 
0 00 00 00 
ae fe(x) COSg(x) tdz, QED. 
0 
(4.10) 
Similarly one can prove the following 
Theorem. The generalized sine (4.8) satisfies the identity 
(sing (2) x)’ =1—fa(zx)sing(,) 2 dz. (4.11) 
0 


Then we prove the validity of 


Theorem. The generalized cosine and generalized sine are linearly 
independent solutions of differential equation (3.1). 


Proof. The Wronskian of two solutions of the differential equation 


(3.1) is . 
W(y1, y2) = 


Y1 Ya 
yi V9 


= Y1Y2 — YY: (4.12) 


From the Liouville theorem of canonical differential equations of the second 
order it is valid that 17(2) = const., which we can choose to be W(2 = 0). 
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On the other hand from the relations (4.5) and (4.7) we have yi(0) = 1, 
y,(0) = 0, yo(0) = 0 and y}(0) = 1. Therefore, from the expression 
(4.12) we obtain that W(0) = 1, so the solutions y, and y2 are linearly 
independent, QED. 
The above theorem leads to the basic relation of generalized trigonom- 
etry of second order 
COSg(x) £ SiNa(x) 2 


(cosa(x) x)’ (sing (2) 2) =i (4.13) 


which in the special case a(x) = 1 transforms into the standard relation 
of the usual trigonometry 
cos* x + sin? z = 1. (4.14) 


Analogously with usual trigonometry one can define generalized tan- 
gent and cotangent 


SINg(x) 2 


tang(2)t = ae (4.15) 
cOta(x) 2 = mee (4.16) 
S11q(2x) x 


which, of course, in the special case a(x) = 1 transforms into the usual 
tangent and cotangent. 

In equation (3.1) let the coefficient be negative, i.e. a(az) = —|a(zx)]. 
Then the series (3.7) and (3.8) converge uniformly and we can analogously 
with the general cosine and sine (4.6) and (4.8) define the general cosine 
hyperbolic and general sine hyperbolic 


yi(2) = cosha(,)% =1 aa ja(x)|dx? + ff |a(x)|dax? ff ja(x)|da? —-.- 
00 00 
(4.17) 


yo(2) =sinhg(z)e = at fzla(x)|dz* + ff |a(z)|dz? ff 2ja(x)|de? —..- 
00 00 00 
(4.18) 


and also generalize tangent and cotangent hyperbolic 
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sinha(z) t 


t = a 

anha(x) £ zoshags) #° (4.19) 
cosh a(x) x 

= ———. 2 

cotha(2) 2 sake (4.20) 
Analogously with the Euler formula 
e” = sinhz+coshz (4.21) 
one can define the generalized exponential function 

€a(x) = Sinha(2) t+ cosha(z) 2 (4.22) 


which the product €7(,. ° Can) # 1, except when e(z) = 1. 
- Finally, we will prove the theorem of the parity of sine and cosine 
function. 


Theorem. For cosine and sine functions the following identity is 
valid : 


COSq(—z)(—2) = COSq(2)(2), sing(—z)(—2) = —singz)(2). (4.23) 


The’Proof is evident according to definition (4.6) 


COSg(—z)(—%) = 1— y dt | a(—t)dt 
ae (4.24) 


+ fdt f a(—t)dt f dt f af—t)dt—--- 
0 0 0 


0 
Then by the replacement —t = q in the right hand side of (4.24) we have 


COSg(—z)(—2) =1— f(—1) dg f a(q)(—1) dg 
+ (-1) aq | a(9)(—1) dg J(—1) da f alq)(-1) agate: 


= 1— {dq fata)dg + f dg falg)dg f dq J a(g)dq —--- 


= 0084(2)(2). a 
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For sine function the proof is similar, see formula (4.8). 

In the special case when the coefficient of the differential equation is 
an even function, i.e. a(—2) = a(z) from the theorem it follows that the 
cosine solution is an even function and the sine solution is an odd function 


COSq(z)(—2) = COSg(2)(Z), sing(z)(—2) = ~ sing(z)(2). (4.26) 


That means the general solution (in a wider sense) of the canonical dif- 
ferential equation of the second order is a sum of odd functions and even 
functions. 

It is evident that both parity relations (4.23) and (4.26) which are 
valid for sine function also are valid for tangent and cotangent functions. 

In in the previous consideration we have introduced the generalized 
trigonometry of canonical differential equation of the second order. But 
these generalized trigonometric functions are formal. They are defined 
only via the multiple integrals (3.7), (3.8) etc. However, when we say 
*trigonometry” we always think of the triangle, cathetus, hypotenuse ... 
The question is: does there exist a connection between rectangle triangle 
and generalized trigonometric function? 

The answer is affirmative [6]. 

First we define what is known as phase” space as a two dimensional 
space where the abscissa and ordinate are a particular solution of the 
differential equation of the second order (3.1) yi(z) and y2(2). The linear 
independence of the solutions lead to the conclusion that the Wronscian 
(4.12) is equal to 1. 

Now instead the orthogonal coordinates y; and y2 we introduce the 
polar coordinates o(2) and y(z) such that 


ys(2) = o(2) cos y(2), 


(4.27) 
y2(2z) = e(x) sin y(z). 
Then from the condition W(y;, yz) = 1 and (4.12) we obtain 
mys — My =1, (4.28) 


or if the equation (4.28) is transformed into the polar coordinates (4.27) 
it is valid that 


dy 
aT 2 
oT, L (4.29) 
After the integration of the expression (4.29) we have 
va) 
| o’(y) dp = 2P, (4.30) 


where P is called the sector surface (Fig. 4.1). So, if one represents the 
particular solutions y; and yo of the canonical differential equation as a 
curve in the phase space, the meaning of the variable z is a double sector 
surface. From the above formula it also follows that dP/dzx = 1/2. The 
sectorial velocity (in the phase space) is a constant! Does it resemble Ke- 
pler’s second law? However Kepler’s law of mechanics is a consequence of 
a motion in three dimensional space in the presence of a conservative cen- 
tral force. Our differential equation corresponds to dimensional harminic 
oscillator with time-dependent mass. It is not clear that there is an equiv- 
alence between these two processes. 


0 ¥ yi (x) 0 1 yi (x) 


Fig. 4.1 The sector surface. Fig. 4.2 The types of phase graphs. 


In the special case when a(x) = 1, the particular solutions (3.7) and 
(3.8) become the usual cosine and sine functions: y; = cosz, y2 = sing. 
Then the phase trajectory is a circle with a radius 9 = 1, and from the 
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integral (4.30) it follows that 2 = f(? dy = y. The meaning of the variable 
x is usual angle. 

Fig. 4.2 shows the global behaviour of the phase trajectory of the 
canonical differential equations of the second order in the neighborhood 
of the point (1,0), for three types: a) elliptic type a(x) > 0, b) hyperbolic 
type a(z) <0 and c) parabolic type a(x) = 0. The tangent on the phase 
trajectory in the point (1,0) is the vertical line. It can be proved simply 
by the calculation of the derivative and using formulae (4.9) and (4.11) 


i SiNg(z) 2)! 
lim La lim (sina(s) 2)” 
t—0 dy, 40 (coSqcz) £)! 


4, 
i 1 — fy a(x) sing(z) & 1 seb (4.31) 
= lim ——.—— "pn — an 
z—0 — fo a(x) cosg(z) £ —0 me 
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| 
| 
| 
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tanha(e) & 


pfs eu 


0 COSg(z) 2 C 1 yi(2) 0 coshg(2)® yi (a) 
Fig. 4.3 Generalized trigonometry. Fig. 4.4 Generalized hyperbolic 
trigonometry. 


Now the geometrical interpretation of the generalized trigonometrical 
function has been shown in the Figs. 4.3 and 4.4. Namely, in Fig. 4.3 in 
the rectangular triangle OBD the sector line OC represents cos,q(2) 2, AC 
represents sin,(,)2 while BD represents tan,(,) . On the other hand, in 
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the hyperbolic case (Fig. 4.4) the line segment OC’ represents the func- 
tion cosh,(,) z, A’C"’ corresponds to sinha.) 2, while B'D' corresponds to 
tanhg(2) &. 

The knowledge of phase trajectory can approach for one kind of in- 
verse problem. Namely, if the form of phase graph @ = e(y) is given, one 
can find [6] the coefficient a(x) and particular solutions of the canonical 
differential equation. 


5. Canonical Equation of n-th Order 


The equation 
y™ +a(z)y =0 (5.1) 


is known as a.canonical differential equation of n-th order. Repeating all 
previotis consideration gives the 
Theorem. The differential equation (5.1) has a particular integral 


va(e) =1— ff fala) de” + f fo fa(w) de” f f+ fale) dz” — +. 
00 0 00 0 00 0 
>_—_ SS SS 


(5.2) 


Proof. After differentiation of the expression (5.2) m times we have 


ys") = -a(2) +a(2) f [--- fale) de” 
00 0 


—a(2) ff fae) de” f f--- fale) de” +... 
00 0 00 0 
—=——— YS — 


=-a(z) |1-ff---fa(x)de” +f f+ fa(n)de” f f+ fale) de® —-. 
00 0 00 0 00 0 
——” —S_—’ 


: 
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1.€. 
yi” +a(z)yi = 0, QED. (5.4) 


By analogy, it can be proved that the functions 


yale) = 2" f fafa" Vale) da +f [++ J ale)de” ff fa" Ya(2) de®—--, 
00 0 00 0 00 0 


(5.7) 
are the particular integrals of differential equation (5.1). The above par- 
ticular integrals can be written in the closed form [5] 


yr(w) = SO(-1)F ff fala) de ---f f+ fa(x) de”, 
k=0 00 0 00 


k n-times integrals 
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(5.8) 


gn(2) = Yt- IT ff Jae) a da” of Foe Fa"Va(e) do”. 
ee —— 


a 


k n-—times integrals 


The characteristics of solutions (5.8) is that in the neighborhood of 
the point z = 0 they have the behavior of a integer degree: yi(z) ~ 1, 
yo(t) ~ 2, ys(t) ~ 2, ++, yn(e) ~ oP}, 

Theorem. The system of the functions (5.8) is a fundamental 
system of particular integrals of canonical equation (5.1). 


Proof. The Wronskian of solutions (5.8) is 


Y1 Y2 tae Yn 
v4 Yo +. Uh 
W (yi, Yars++s Yn) = : ; ‘ : 2 (5.9) 


yr) yr) ate yr 


If we denote Wo = W[y1(0), y2(0),.--, ¥n(0)] then by an elementary 
approach to the formula (5.8) there follows the result 


1 0 O 0 0 
0 1! O 0 0 
w,=|9 O 2! 0 0 =(n—1)!---2!1!40 QED. 


0°20: Oo gx 0 Gay 
(5.10) 
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So, the general solution y(x) of the canonical differential equation of 
n-th order is given by the formula 


y(z) = Cry + Cayo +--+: + Cnyn, (5.11) 


where Cj, 1 = 1,2,...,n, are arbitrary constants, while y; are the funda- 
mental solutions (5.8). 

Theorem. The series (5.8) are convergent uniformly for every 
analytical coefficient a(z). 


Proof. From |a(z)| < M we have that 


, es deala = 
1 Le i eee en eaey 
at ee », aby 


1 (|x|) 


je|"** 2 lef" - 
< M M* 7 rr = = eho 
Iya] < |x| + (n +1)! ot (2n +1)! a 2 YM (nk+1)! 
Bier 


|x|?" 


(2n — 1)! (in 1)! 


ie (vie) 
oe ey my oe 


lyn| < el" + M 


k=1 

(5.12) 
and because the majoring series is convergent (the D’Alamberts test), the 

series (5.8) converge uniformly, QED. 
Similarly to previous considerations, we can study the phase space 
of the canonical differential equation of n-th order (5.1). It is an n- di- 
mensional Euclidean space of the particular solutions (5.8). Assuming the 
constants Cy, C2, ... Cy in (5.11) as scale factors one can represent the 
solution of differential equation (5.1) by the point in this n dimensional 
space. Then the evolution of the solution y(z) is a curve in this space. 
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The elementary n-volume of the n-body constructed on the curve 7(2) in 
n dimensional phase space is 


1 
dv™ = 5 Vly (x), yo(2),.--,Yn(x)] de. (5.13) 


where W is the Wronskian (5.9). Then from the relations (5.10) and 
(5.13) we obtain the geometrical interpretation of the parameter x of the 
canonical differential equation of n-th order [6] 


(n) 
ee (5.14) 
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